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Abstract 

The performance of two recent ab initio computational thermochemistry 
schemes, Wl and W2 theory [J. M.L. Martin and G. de Oliveira, J. Chem. 
Phys. Ill, 1843 (1999)], is assessed for an enlarged sample of thermochemical 
data consisting of the ionization potentials and electron affinities in the G2-1 
and G2-2 sets, as well as the heats of formation in the G2-1 and a subset of the 
G2-2 set. We find Wl theory to be several times more accurate for ionization 
potentials and electron affinities than commonly used (and less expensive) 
computational thermochemistry schemes such as G2, G3, and CBS-QB3: W2 
theory represents a slight improvement for electron affinities but no signifi- 
cant one for ionization potentials. The use of a two-point A + B /L 5 rather 
than a three-point A + B /C L extrapolation for the SCF component greatly 
enhances the numerical stability of the Wl method for systems with slow 
basis set convergence. Inclusion of first-order spin-orbit coupling is essential 
for accurate ionization potentials and electron affinities involving degenerate 
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electronic states: inner-shell correlation is somewhat more important for ion- 
ization potentials than for electron affinities, while scalar relativistic effects are 
required for the highest accuracy. The mean deviation from experiment for 
the G2-1 heats of formation is within the average experimental uncertainty. 
Wl theory appears to be a valuable tool for obtaining benchmark quality 
proton affinities. 
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I. INTRODUCTION 



Development of models based on molecular orbital theory for theoretical thermochemistry 
involves five key steps : defining a target accuracy, formulation of theory, implementation 
through programs, validating the models against reliable experimental values and predic- 
tion on any molecular system by the end user. At present, only ab initio methods can 
claim "chemical accuracy" (commonly defined as 1 kcal/mol) for small and medium sized 
molecules. The most popular such methods are the Gaussian- n (Gn) theories [@-0 of Pople 
and coworkers (which are based on a combination of additivity approximations and empirical 
corrections applied to relatively low-level calculations), followed by the CBS approaches |5|-[7| 
of Petersson and coworkers which are intricate combinations of extrapolation and empirical 
correction schemes. 

Very recently, Martin and de Oliveira presented two theoretical thermochemistry schemes 
known as Wl and W2 (Weizmann-1 and Weizmann-2) theory ||, which aim at 'benchmark 
accuracy', defined by these authors as 1 kJ/mol (0.24 kcal/mol). For a set of 28 experimen- 
tally very precisely known molecular total atomization energies, the more cost-effective of 
the two schemes, Wl theory, achieved a mean absolute error of 0.37 kcal/mol, while the more 
rigorous of the two schemes, W2 theory, achieved a mean absolute error of 0.23 kcal/mol. (It 
should be pointed out that these methods are free of parameters derived from experiment: 
Wl theory does contain one parameter — the exponent for the valence correlation extrap- 
olation — that is derived from W2 calculations, not experiment.) Martin later proposed 
a minor modification of Wl theory, denoted Wl' theory fH, which appeared to yield con- 
siderably improved results for second-row compounds at no additional cost. (For first-row 
compounds, it is identical to Wl theory.) 

In recent years, density functional theory (DFT) methods have also reached a level of 
sophistication where they can provide thermochemical data to within a few kcal/mol, notably 
the hybrid B3LYP (Becke 3-parameter exchange with Lee- Yang-Parr correlation [|il]| ,|Tl"|j ) 
and B97 (Becke 1997 |L2]]) exchange-correlation functionals, but also the 'pure DFT' HCTH 
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(Hamprecht-Cohen-Tozer-Handy fl3| ) exchange-correlation functional. A recent collection 
of reviews on computational thermochemistry methods may be found in an ACS Symposia 
volume edited by Irikura and Frurip. [JT4"| 

The most fundamental thermochemical property of a compound, from an experimental 
point of view, is its heat of formation (AHj- ) in the gas phase. From a computational chem- 



istry point of view, the total atomization energy (TAE, TjDo) |15 is the most fundamental 
such quantity. Using the experimental heats of formation of the atoms in the gas phase, 
TAEs can be directly related to the gas-phase heats of formation. 

Prior to proper application of any new model by the end user, it should be tested against 
known high quality experimental results. For this purpose, Pople and coworkers proposed 
two standard test sets of thermochemical data: the G2-1 test set being the smaller and 



containing small molecules, and the G2-2 test set |[16|Ji7[| containing larger systems. These 
sets of thermochemical data, covering 148 neutral and 146 ionic species, have been used 
fairly extensively (e.g. [18|]) to test the performance of various computational thermochem- 
istry methods, notably the Gn theories and their variants [ |TB] , P~7| , [TD^ P^] , density functional 
methods [gjrj^g , and the CBS family of methods 

The main problem with the G2-1 and G2-2 test sets for heats of formation is the limited 
accuracy of the experimental data themselves. These were critically reviewed by Liebman 
and Johnson |26| , who concluded that less than half of the data even met the less rigorous 1 
kcal/mol accuracy criterion. For methods of the W1/W1'/W2 type, this is clearly a major 
impediment to their validation for a larger experimental data set, and alternatives need to 
be sought. 

Pople and coworkers also defined G2-1 || and G2-2 |17j data sets for ionization potentials 
and electron affinities. The accuracy of these experimental data is much more satisfactory, 
and it could be argued that they are in fact more suitable test sets for the validation of 
theoretical thermochemistry methods. While ionization potentials are comparatively easy to 
reproduce well, electron affinities are a very taxing test for any electronic structure method. 
This is true both in terms of the basis set (addition of the electron entails a profound change 
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in the spatial extent of the wave function) and in terms of the electron correlation method 
(effectively, the number of interacting electrons in the system changes). 

The purpose of the present paper is to assess the performance of Wl and W2 theory for 
an extended data set of thermochemical data. Specifically, for Wl theory we shall consider 
the G2-1 and G2-2 datasets for ionization potentials and electron affinities, as well as the 
G2-1 and a subset of the G2-2 data set for heats of formation. W2 theory will be considered 
for the G2-1 data set for IPs, EAs, and heats of formation. Finally we shall turn to Wl and 
W2 theory for proton affinities. 



II. COMPUTATIONAL METHODS 

Geometry optimizations and vibrational frequency calculations using the B3LYP (Becke 
3-parameter-Lee- Yang-Parr JT^JTl|) density functional method have been carried out using 
Gaussian 98 revision A7 |27j . (Following the recommendations in Ref . |28| larger grids than 
the default were used in the DFT calculations if necessary, specifically a pruned (99,590) 
grid for integration and gradients, and a pruned (50,194) grid for the solution of the coupled 
perturbed Kohn-Sham equations.) All other calculations were carried out using MOLPRO 



98.1 |2{J and a driver for the W1/W2 calculations |30|] written in MOLPRO's scripting 
language. The lion's share of the calculations was carried out on a Compaq ES40 with four 
667 MHz Alpha EV67 CPUs, and a scratch volume consisting of six 18GB SCSI 2 ultrawide 
disks striped in software. Remaining calculations were carried out on the SGI Origin 2000 
of the Faculty of Chemistry. 

The SCF and valence correlation basis sets are Dunning's augmented correlation consis- 



tent n-tuple zeta [pTH 33|| (aug-cc-pVnZ) basis sets; for second-row atoms, high-exponent d 



and / functions were added (denoted '+2d' or '+2dlf ') as was found repeatedly p4] , |35| to be 
necessary for accommodating inner-shell correlation effects. (Unless indicated otherwise, reg- 
ular 'un-augmented' cc-pVnZ basis sets are used for H, Li, and Be, as well as cc-pVnZ+2dlf 
for Na and Mg.) For all remaining steps (inner-shell correlation, scalar relativistic effects, 
and spin-orbit coupling), the MTsmall (Martin- Taylor small M) core correlation basis set 



was used. Restricted open-shell wave functions were used throughout for open-shell species. 

The W1/W2 energy consists of seven components, each of which we shall detail in turn 
for reasons of clarity and self-containedness. 

(0) Reference geometries are obtained at the B3LYP/cc-pVTZ+l level in the case of Wl 
theory and at the CCSD(T)/cc-pVQZ+l level in the case of W2 theory. In both cases, 
the '+1' signifies the addition |36| of a high-exponent d function to second- row elements, 



the exponent having been set equal to the highest d exponent in the corresponding cc-pV5Z 
basis set. 

(1) In the original Wl and W2 papers, the SCF limit was obtained by geometric extrapo- 



lation 37 



TAE SCF (n) = TAE SC F,oo + A.B n (1) 

of the molecular total atomization energy TAE computed using cc-pVnZ+2dlf basis sets, 
where for W2 theory n = {T, Q, 5} (with Z={3,4,5}) , and for Wl theory n = {D,T,Q} 
(with Z={2,3,4}). (In practice, this means that TAEscf,oo = TAEscf,u — (TAEscf,u — 
TAE SC F,n-i) 2 /(TAE S cF,n ~ 2TAE SCF ,n-i + TAE SCF ,n~2)- The exponential convergence 
behavior of the SCF energy has repeatedly been demonstrated empirically (e.g. by Jensen 
p8| and by Martin and Taylor ||39|| ) in comparisons with numerical Hartree-Fock energies 
for small molecules. 

The geometric formula, for this particular application, has the minor disadvantage that 
its extrapolated limit depends on whether the extrapolation is carried out on TAE, or on the 
constituent energies. (In practice the differences are quite minor.) Based on the asymptotic 
convergence behavior ||40|— 42 of the pair energy in an electron pair that does not have an 



interelectronic cusp, Petersson and coworkers |E| previously considered (within the context 
of their CBS family of methods) an alternative formula E SCF (n) = E S cf,oo + YliZ n +i • 
(Using Euler-Maclaurin summation, we find this to be equivalent to the simple two-point 
formula Sscf(^) = Escf,oo + A/n 5 + 0(n~ 7 ).) Martin and Taylor |44| previously considered 
the difference between the three-point geometric formula (with n={T,Q,5}) and the two- 
point formula (with n={Q,5}) for a small set of molecules and found the differences to be 
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negligible. We also find this to be the case for the much larger sample of molecules surveyed 
here. 

However, when considering three-point geometric n={D,T,Q} versus two-point n={T,Q} 
in the present work, we found that, while the differences are quite small for almost all first- 
row and most second-row systems, significant differences (in excess of 1 kcal/mol) exist 
for a few first-row systems (e.g. LiF) and a rather larger number of second-row systems 
(e.g. many silicon compounds). Some conspicuous examples can be found in Table I. Upon 
closer inspection, this was revealed to be caused by the three points lying nearly on a 
straight line, causing the geometric extrapolation to yield erratic results. The two-point 
extrapolation is invariably closer to the extrapolated limit obtained from the larger basis sets: 
in unproblematic cases, it yields essentially the same results as the three-point extrapolation. 
As a result, we are recommending that the two-point A + B/l 5 extrapolation be used in Wl 
and W2 theory from now on: one beneficial side effect is that the extrapolated limit for this 
two-point formula is easily seen to be independent of whether the extrapolation is carried 
out on the molecule or the constituent atoms. Since the SCF component was the only 
component for which such an ambiguity existed in the original Wl and W2 theory, this 
permits the quoting of "total Wl and W2 energies" for arbitrary systems. 

It is also seen in Table I that the main argument in favor of Wl' theory (in which the 
AVTZ+2dlf basis set is replaced by an AVTZ+2d basis set, for balance reasons ||) over 
standard Wl theory, namely an SCF limit in better agreement with that obtained from 
larger basis sets, appears to be obviated by the new extrapolation. We shall not consider 
Wl' theory further in the course of this paper. 



(2) In the W2 case, the CCSD (coupled cluster with all singles and doubles (45 



va- 



lence correlation contribution to TAE is obtained using the aug-cc-pVQZ+2dlf and aug- 
cc-pV5Z+2dlf basis sets, then extrapolated to the infinite basis limit using the expression 
g§ E(l) = E 00 + A/L 3 . (In practice, this means = E t + (E t - £,_i)/((^) 3 - 1).) 
The arguments in favor of this expression (derived from the known asymptotic conver- 
gence behavior of the interelectronic cusp [j40H42|1) have been detailed at length elsewhere 
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|8| J46| . |47[| and will not be repeated here. In the Wl case, the unmodified expression leads 
to systematically overestimated correlation contributions to TAE ||: here we employ 
= Ei + (Ei - ^_i)/(( I ^ T ) a - 1), where a = 3.22 was determined § to yield the 
best agreement with the extrapolated W2 CCSD energies. Both for Wl and W2 theory, the 
largest basis set CCSD calculation is carried out (except for very small systems) using the 
direct CCSD algorithm of Schutz, Lindh, and Werner f48| as implemented in MOLPR098.1. 

(3) The contribution of connected triple excitations is obtained at the CCSD(T) level (CCSD 
with a quasiperturbative a posteriori correction for connected triple excitations fl49fl). As the 
T3 contribution is known to converge more rapidly than the contribution of exp(7\ + T2), 
hence this contribution is obtained from CCSD(T) calculations with the smaller two basis 
sets and extrapolated to the infinite-basis limit using E^ = E[ + {E\ — / ((rzr) a — 1), 
where, as for the CCSD energy, a=3 exactly for W2 theory and a=3.22 for Wl theory. (For 
open-shell systems, the definitions of the restricted CCSD and CCSD(T) energy as given in 
Ref. [|5T] has been used.) 

(4) The inner-shell correlation contribution is computed as the difference between 
CCSD(T)/MTsmall |J values with and without constraining the inner-shell orbitals to be 
doubly occupied. (In the case of the second-row elements, the very deep-lying (ls)-like 
orbitals are constrained to be doubly occupied throughout.) 

(5) The scalar relativistic contribution is computed as expectation values of the one-electron 
Darwin and mass- velocity (DMV) operators []52| , |53| for the ACPF/MTsmall (averaged cou- 
pled pair functional |5^| ) wave function, with all inner-shell electrons correlated except the 
(ls)-like orbitals of second- row elements. Bauschlicher |55 demonstrated that, for first- 
and second-row systems, this approach yields essentially identical results to more rigorous 
relativistic calculations. 

(6) For closed-shell systems, or open-shell systems in nondegenerate electronic states, there 
is no molecular first-order spin-orbit contribution, and the contribution to TAE is merely 
the sum of the atomic spin-orbit corrections. For open-shell systems in degenerate states, 
we have calculated spin-orbit corrections at the CISD level with the MTsmall basis set, and 
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again correlating all electrons except for the (Is) on second-row elements. 
(7) The molecular zero-point energy and thermal corrections were obtained at the 
B3LYP/cc-pVTZ+l level. The zero-point energies within the harmonic approximation are 
scaled by 0.985, primarily to correct for anharmonicity. The scale factor was obtained || by 
comparison with experimental (or high-level theoretical) anharmonic zero-point energies for 
28 molecules. 

(Adiabatic) electron affinities (EAs) are calculated as the difference between the TAEo 
values of the anion and the corresponding neutral species, at their respective optimized 
geometries: 

EA = TAE (anion) - TAE (neutral) (2) 

Likewise, the (adiabatic) ionization potentials (IPs) are calculated as the difference in total 
atomization energies at K of the cation and the corresponding neutral, at their respective 
optimized geometries: 

IP Q = TAE (neutral) - TAE (cation) (3) 

Theoretical heats of formation at K were calculated by subtracting the Wn calculated 
TAE (£T>o) value from experimental enthalpies of formation of the isolated atoms. For any 
molecule, such as A^B^H^, the heat of formation at K is given by 

AH° f (A x B y H z , OK) = xAH^A, OK) + yAH° f (B, OK) + zAH° f (H, OK) - ED Q . (4) 

The CODATA |56| values of the atomic AH® are used with the exception of boron and 



silicon, for which we have used revised values recommended by Bauschlicher, Martin, and 



Taylor |57|j58|1 for boron and by Martin and Taylor |59[ for silicon. Theoretical heats of 



formation at 298 K (AHj) are calculated by correction to AH® as follows: 

AH° f (A x B y H z ,298K) = AH° f (A x B y H z ,0K) 

+ [H°(A x B y H z , 298K) - H°(A x B y H z , OK)} 

- xhcf 298 [A, st] - yhcf 2 9 8 [B, st] - zhcf 298 [H, st] (5) 



While the enthalpy functions he fx = Ht — Hq for the molecule are obtained using the RRHO 
(rigid rotor-harmonic oscillator) approximation from the unsealed B3LYP/cc-pVTZ+l vi- 
brational frequencies, the enthalpy functions for the standard states of the elements are taken 
directly from CODATA. 

Proton affinities (PA) are obtained from the total atomization energies at K as follows: 

PA (B) = TAE (BB. + ) - TAE (B) (6) 

Finally, PAs at 298 K are calculated by correction to PAo as follows: 

PA 298 (B) = PA (B) + hcf 298 (M+) - hcf 298 (B) - — (7) 

where 5RT/2 is the enthalpy function of the H + ion. 



III. TEST SETS USED 

The original G2-1 ion test set consists of 25 EAs and 38 IPs while the G2-2 test set 
included 33 EAs and 50 IPs. In the G3 paper, Curtiss et al. |4[] applied G3 theory to the 
G2-1 and G2-2 test sets, minus three ionization potentials (due to the size of the molecules 
concerned). In this study we exclude five additional ionization potentials and one electron 
affinity from the G2-2 test set, for the same reason. Both Wl and W2 theories were evaluated 
for the G2-1 test set, while only Wl theory was considered for the G2-2 test set. For the 
purpose of evaluation of Aif^, the original G2-1 neutral test set consists of 55 molecules 
while the G2-2 test set includes 93 molecules. Again, we have considered both Wl and W2 
theories for the G2-1 neutral test set of molecules and compared the results with G2, G3 and 
CBS-Q values. (For the W2 calculations, three species were omitted because of their size.) 
It should be noted that the G2-2 test set contains several fairly large molecules and some 
of the experimental AHj for the species in G2-2 test set possess large uncertainties as well 
as several experimental values spanning a wide range. Therefore, we have selected a subset 
of 27 out of the 93 G2-2 neutral molecules, which are tractably small and for which the 
experimental enthalpies are reasonably accurate. To these molecules we applied Wl theory, 
and to a subset of them W2 theory. 
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IV. RESULTS AND DISCUSSION 



A. Electron Affinities 

We shall first consider the G2-1 test set. A breakdown of the different components of the 
Wl values is given in Table II, while a comparison between various levels of theory (including 
Wl, W2, G2, G3, and CBS-Q) and experiment is given in Table III. 

We note first that a substantial number of species have negative electron affinities at 
the SCF level: the binding of the electron results from the additional correlation energy in 
those cases. Inclusion of connected triple excitations is essential. In contrast, inner-shell 
correlation does not appear to be very important for the G2-1 EAs. Scalar relativistic effects 
are somewhat more important: with the exception of CH 2 , CH 3 , and SiH 3 they uniformly 
decrease the electron affinity. As expected, the scalar relativistic effect is somewhat more 
important in second-row than in first-row systems. The change in the zero-point energies 
can be fairly substantial, particularly for hydrides. 

For heats of formation of closed-shell systems (and open-shell systems with nondegenerate 
ground states), the molecular first-order spin-orbit splitting vanishes, reducing the spin-orbit 
correction to the sum of the corrections for the constituent atoms. Since a fair number of the 
species in the G2-1 set have a degenerate state for either the neutral or the anionic system (or 
in fact for both), some account for molecular spin-orbit coupling cannot be avoided. We have 
considered a number of (relatively) inexpensive approximations within the MTsmall basis 
set used for the core correlation and scalar relativistic steps, including SCF, CISD, and CISD 
with inner-shell electrons correlated. The computed corrections at these levels of theory for 
a number of (neutral, cationic, and anionic) species have been compared in Table IV with 
values obtained from experimental fine structures. For most first-row species, satisfactory 
results are already obtained at the SCF and definitely at the CISD level; for the second- 
row species, correlation from the (2s2p) inner-shell orbitals appears to be essential, as was 
previously found by de Oliveira et al. |6(| for the 2nd-row atoms and by Nicklass et al. |Tl| 
in a convergence study for the halogen atoms. (The rather weak basis set dependence found 
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by these latter authors |61| is consistent with our own findings.) Only for the CIO molecule 
we find a substantial error: inspection of the spectroscopic constants for the few lowest 
states reveals that the X 2 I1 ground state in fact has a anomalously small splitting compared 
to the A 2 I1 state; since this latter state mixes in quite prominently into the X 2 I1 wave 
function, the splitting is severely underestimated unless the A 2 I1 state is admitted to the 
zero-order wave function. Using a CASSCF reference space consisting of the valence orbitals 
except for the CI (3s) and 0(2s) like orbitals, and supplemented with the first Rydberg tc 
orbitals, yields a spin-orbit correction in excellent agreement with experiment. 

As can be seen from Table II, these spin-orbit corrections are in fact essential for good 
agreement with experiment for several of the systems. At the Wl level, we find a mean abso- 
lute discrepancy (MAD) from experiment of 0.016 eV (Table III), which is a quite substantial 
improvement over the G2 (0.057 eV), G3 (0.049 eV), and CBS-QB3 (0.054 eV) values. Per- 
haps even more importantly, the maximum error is likewise much smaller, 0.051 eV for CH3, 
followed by 0.043 eV for SiH 2 . In the case of CH 3 , not only is the electron affinity very 
small (G3 and CBS-QB3 in fact predict the wrong sign), but the harmonic approximation 
for the zero-point energy is of dubious reliability (see Schwenke |63[| ). W2 represents a minor 
improvement over Wl, at vastly greater computational expense: MAD=0.012 eV. Using 
even larger basis sets, de Oliveira et al. found the mean absolute error for the atoms H, 
B-F, and Al-Cl to be 0.009 eV at the CCSD(T) level; by employing CCSDT and full CI 
corrections, this error could be reduced by an order of magnitude. (The importance of these 
corrections was about evenly split between higher-order T 3 effects and effects of connected 
quadruple excitations, T 4 .) We conclude that the accuracy of W2 theory (and, to a lesser 
extent, Wl theory) is mostly determined by the imperfections in the CCSD(T) method. 

We shall now consider the G2-2 set of electron affinities, for which only Wl (not W2) 
calculations were carried out. A comparison with other theoretical thermochemistries and 
with experiment is given in Table V, while a breakdown of contributions is given in Table 
S-I [0. 

The trends seen for the G2-1 set largely continue for the G2-2 set. However, inner-shell 
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correlation is somewhat more important for some species (e.g. Al, because of the small 
subvalence/valence gap, C 2 , and S 2 0). One exception to the general trends is that electron 
correlation in fact decreases EA(C2): this is an artifact of the multireference character of 
the X X E+ state. At first sight, scalar relativistic effects seem less important, but this is an 
artifact of the relative preponderance of first-row species compared to the G2-2 set. 

Standard Wl results for Li and Na would not involve diffuse functions on these low- 
electronegativity atoms. Not surprisingly, very poor electron affinities are thus obtained. 
We have optimized diffuse functions (available in the supplementary material) for Li, Be, 
Na, and Mg to accompany the standard cc-pVnZ basis set: the exponents were optimized 
individually for each angular momentum at the CISD level for the atomic anion. The Wlaug 
results obtained with these basis sets are in near-perfect agreement with experiment. (This 
can reasonably be expected since the electron correlation methods used are exact within the 
finite basis set for the valence correlation contributions.) 

As for the G2-1 set, we find Wl theory to be quite substantially more reliable than G2 
and G3 theory. Substantial discrepancies between Wl theory and experiment are found for 
ozone, CH 2 NC, and FO: the first two species (and, to a lesser extent, FO) exhibit strong 
nondynamical correlation effects, and hence methods that do not include corrections for cor- 
relation effects beyond CCSD(T) are expected to yield poor results. G3 theory fortuitously 
agrees better with experiment than Wl theory for these species. 

B. Ionization Potentials 

We shall again first consider the G2-1 set. A breakdown of components in the Wl 
computed values can be found in Table VI, while a comparison with experiment and less 
expensive theoretical thermochemistry methods can be found in Table VII. The relative 
importance of correlation is smaller than for the electron affinities: yet in absolute terms its 
contribution is almost as significant as in the EA case. While connected triple excitations 
appear to be somewhat less important than for EAs, they can certainly not be neglected 
with impunity. Inner-shell correlation contributions, on the other hand, are more important 
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than in the EA case because the valence excitation creates a "hole" into which core electrons 
can be excited. The large contributions for Na (4 kcal/mol) and Mg (2 kcal/mol) come as 
no surprise given the small core-valence gap in these atoms. Scalar relativistic contributions 
are important for accurate work: with the exception of Li-C and Na-Mg, they consistently 
lower the IP. Like for the EAs, we see substantial zero-point effects for the hydrides: in the 
case of CH4, this contribution is especially large because of the known fluxional nature |>l| of 
the CH^ cation. And again, spin-orbit splitting is a factor to be reckoned with, particularly 
for such second-row species as exhibit first-order spin-orbit splitting. 

Agreement with experiment is highly satisfactory at the Wl level, except for CH4 where 
an atypically large discrepancy is seen. Upon inspection, it is revealed that the B3LYP/cc- 
pVTZ geometry for CH4 is qualitatively incorrect, exhibiting D 2 rather than C 2v symmetry. 
This reflects itself both in an error in the total energy for the cation and in an error in the 
zero-point contribution. Using a CCSD(T)/cc-pVTZ reference geometry, excellent agreement 
with experiment is in fact obtained. 

In an attempt to ascertain whether this issue is specific to the B3LYP exchange- 
correlation functional, we carried out geometry optimizations and vibrational frequency 
calculations for CH4 using the cc-pVTZ basis set and a variety of exchange-correlation 
functionals, including B3P86 (Becke 3-parameter exchange with Perdew-1986 correlation 



), B3PW91 (Becke 3-parameter exchange with Perdew-Wang-1991 correlation |67]), 
mPWlPW91 (modified Perdew-Wang |)!|), mPWIK (Truhlar's empirical modification of 
the latter @), BHLYP (Becke half-and-half exchange [flj with LYP correlation), BHPW91 



(ditto with PW91 correlation), and BLYP (Becke 1988 exchange |71] with LYP correlation). 
Results are summarized in Table VIII. Only the functionals with 50% Hartree-Fock exchange 
(BHLYP, BHPW91) or nearly so (mPWIK) find this structure to be a local minimum, while 
all other functionals find an imaginary frequency of a 2 symmetry. Following the latter down- 
hill leads to the D 2 d structure. Given that this behavior persists with a fairly wide variety 
of correlation functionals, the problem appears to reside in the exchange functional. (Note 
that since Wl theory does not contain any parameters that depend on the level of theory 
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for the reference geometry, it can quite well be carried out from, say, a mPWIK/cc-pVTZ 
reference geometry for systems which exhibit this type of problem.) 

Mean absolute deviation for Wl is a factor of 3-4 smaller than for the inexpensive 
methods. In this case, only marginal improvement is seen upon going to the much more 
expensive W2 method, which is easily understood in terms of the faster basis set convergence 
for the cation compared to the anion. Again, we have reason to believe that the principal 
factor limiting the accuracy of our calculations are small deficiencies in the CCSD(T) electron 
correlation method. 

Let us now consider the G2-2 ionization potentials (Tables IX and S-II). Most systems in 
that set do not exhibit first-order spin-orbit splitting, the main exceptions being Ne, Ar, OCS, 
and CS2. Most trends from the G2-1 set are continued: one clear exception to the general rule 
is CN, for which electron correlation reduces the IP. At first sight, this system also exhibits a 
large discrepancy of 0.27 eV between theory and experiment, and discrepancies for the more 
approximate G2 and G3 methods are similarly high. An explanation in terms of the extreme 
multireference character of the CN + cation would be tempting: however, we repeated the W2 
calculation using full valence CAS-ACPF instead of CCSD(T) at every step, and found an 
increase by 0.04 eV in our computed value (to 13.93 eV). Upon closer inspection, it appears 
that the 'experimental' IP(CN)=13.6 eV is in fact a propagated transcription error from Ref. 
f72fl . The only experimental value without a large error bar, 14.03±0.02 eV, was obtained 
by Berkowitz et al. f73| from photoionization data for HCN-^H++CN+e" (19.00±0.01 eV) 
and HCN^H+CN++e- (19.43±0.01 eV), as well as the well-established IP(H). The error 
bars on the extrapolated ionization limits may be somewhat optimistic: in addition, it is well 
known from benchmark calculations (e.g. ||74j| : see also [fT5|) that CN - has a very low-lying 
a 3 I1 state (T e =880±100 cm -1 |T4]], or 0.11 eV). It is not inconceivable that the Berkowitz 
et al. value in fact corresponds to generation of the 3 II state (especially since the ground 
state of HCN + is X 2 II): our Wl and W2 computed IPs in that case are both 14.07 eV, in 
excellent agreement with the Berkowitz value. This problem merits further investigation. 

For the purpose of assessment of the error statistics of the various methods, however, 
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we have removed IP(CN) from the sample. Significant discrepancies — out of character 
with the other results — are then still seen for B 2 H 4 , sec-C 3 H 7 , Si 2 H 6 , and CH 3 OF. None of 
these species exhibits severe nondynamical correlation, and we note that there are significant 
discrepancies between G3 and experiment for all these species except Si2H 6 . We suggest that 
these experimental values may need to be reconsidered. In the case of B 2 H 4 , it has previously 
been suggested JF|| that the surface between the D 2 d structure and the doubly bridged C 2v 
isomer is flat enough that rigid molecule treatments may not be appropriate. The remaining 
three cases contain internal rotations, likewise casting doubt on the applicability of the 
RRHO approximation for the zero-point energy. (For both B2H4 and Si2Hg, W2 calculations 
were feasible, and were found to yield essentially the same result as the Wl calculation.) 
Neither neutral nor cationic systems exhibit appreciable multireference character which could 
negatively affect the quality of the Wl and W2 results. Upon eliminating the four doubtful 
species, we find a MAD for the G2-2 test set that is only slightly higher than for the G2-1 
set. Regardless of whether these four species are eliminated, it is clear that Wl represents a 
significant improvement over G2 and G3 theory. 



C. Heats of Formation 

We shall finally turn to heats of formation for a larger set of molecules than was consid- 
ered in the original W1/W2 paper. A comparison with experiment and more approximate 
methods can be found in Table X, while a breakdown by components of the atomization 
energies of the G2-1 set of neutral molecules is given in Table S-III. 

First of all, we note that the mean uncertainty for the experimental values is itself 0.6 
kcal/mol. In fact, the MAD values for Wl and W2 theory stand at 0.6 and 0.5 kcal/mol, 
respectively, suggesting that these theoretical values are in the same reliability range as the 
experimental data. The MAD for W2 theory is more than twice as large as in the original 
W1/W2 paper: however, comparisons there were made against a smaller set of molecules for 
which the experimental uncertainties were all 0.25 kcal/mol or less, mostly 0.1 kcal/mol or 
less. 
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For ten species do discrepancies between Wl theory and experiment reach or exceed 1 
kcal/mol. Out of these, four experimental values carry uncertainties of 1 kcal/mol or more, 
and can be ignored. Of the remaining six, the experimental heat of formation of SiH4 contains 
an ambiguity |f77j , P2 is a notoriously difficult molecule |7§| and carries an uncertainty of 
0.5 kcal/mol, and CIO is strongly multireference and carries an uncertainty of 0.5 kcal/mol. 
(For this latter molecule, however, 'upgrading' the calculation to W2 theory reduces the 
discrepancy with experiment to 0.5 kcal/mol, suggesting that slow basis set convergence 
may be at stake here.) As for Si2Hg, Wl, W2, G2, and CBS-Q theories exhibit similar 
discrepancies from experiment (G3 a somewhat smaller one), strongly suggesting that the 
experimental value may be in error. 

We note that W1/W2 and the less expensive methods "err" in the same direction for P 2 
and CIO as well, suggesting that some revision of the experimental data may be in order there 
as well. Revisions for BeH and NH2 were suggested previously ||79|| : all methods unanimously 
suggest the PH 2 value to be in error. 

As pointed out repeatedly fl50"| , the JANAF heat of formation for SiH 4 is in fact the 
older Gunn and Green value ]5T|] increased by a somewhat arbitrary term of 1 kcal/mol for 
a Si(amorphous)— >Si(cr) phase change. The Wl and W2 results, like a previous benchmark 
study [j32|l , favor the older Gunn-Green value. 

Finally, we selected 26 species with relatively small error bars out of the 93 molecules 
in the G2-2 test set for heats of formation. The experimental Aif J along with deviation of 
Wl, G2, and G3 values from experiment are presented in Table XI, while a summary of our 
computed TAEs and their different components for the 27 G2-2 neutral test molecules set 
calculated at the Wl level is presented in Table S-IV. 

The average discrepancy between the Wl values and experiment for this subset is 0.7 
kcal/mol, compared to an average experimental uncertainty of 0.4 kcal/mol. In order to 
establish the reason for some of the discrepancies, we have carried out W2 calculations on 
selected systems. 

In the case of BF 3 and CF 4 , the culprit appears to be slow basis set convergence in 
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these highly ionic systems. We were unable to complete the CCSD/aug-cc-pV5Z calculation 
for CF 4 : but applying the Wl and W2 extrapolations to the published CCSD(T)/aug-cc- 
pVnZ (n=3,4 and 4,5, respectively) total energies for CF4 of Dixon et al. |83|, we found the 
estimated W2 TAE to be 1.5 kcal/mol lower than the estimated Wl TAE. This accounts 
for essentially all the discrepancy between experiment and Wl theory for CF 4 . (A similar 
phenomenon was previously noted for BF3 ||58|| .) The NO2 molecule exhibits strong nondy- 
namical correlation effects, and the W2 result is actually further removed from experiment 
than the Wl result. C1NO likewise exhibits substantial nondynamical correlation, and here 
the W2 result is basically identical to its Wl counterpart. Improving agreement with experi- 
ment for these two molecules will certainly require accounting for correlation effects neglected 
at the CCSD(T) level. While this may be true to a lesser extent of the N 2 molecule, an 
error in the experimental value cannot completely be ruled out there. 

For F2O, the discrepancies of -0.8 (Wl) and -1.0 (W2) kcal/mol with experiment, as well 
as the absence of significant nondynamical correlation, suggest that the experimental value 
may need to be reinvestigated. 

At least for some of the larger systems (as well as those which have internal rotations), 
neglect of anharmonicity in the zero-point energy may account for part of the discrepancy 
with experiment. 



D. Proton affinities 

Curtiss et al. @], in the original G3 paper, considered proton affinities of eight molecules 
as well. We have calculated Wl and W2 proton affinities for the same systems. However, 
rather than the somewhat older experimental data used by these authors, we have taken 



the data from the very recent compilation by Hunter and Lias |84|. Various methods are 
compared with experiment in Table XII, while a breakdown of the different components at 
the Wl level is given in Table S-V. 

Since protonation/deprotonation is an isogyric reaction, the proton affinity converges 
considerably more rapidly with the level of theory than, say, the heat of formation. (This 
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is expressed, for instance, in the comparatively small contribution of valence correlation, 
and the quite small contributions of inner-shell correlation and scalar relativistic effects.) 
Indeed, we note that only minute differences exist between the Wl and W2 proton affinities: 
Wl theory can basically be considered converged for this purpose. Mean absolute deviation 
from experiment is 0.43 kcal/mol (compared to 1.2 kcal/mol for the inexpensive G3 theory). 
While only a few of the values in Hunter and Lias carry explicit error bars (e.g. water, ±0.7 
kcal/mol, H 2 S, ±1.3 kcal/mol), it is clear that the uncertainty on the computed Wl values 
is considerably lower than that of the experimental values themselves (with the exception of 
H 2 , for which the Hunter and Lias value is a theoretical one). As such, Wl theory should be 
a powerful tool for obtaining benchmark-quality proton affinities: for application to larger 
systems, the inner-shell correlation and scalar relativistic steps can fairly safely be skipped 
for this application. 

The somewhat surprising difference of 0.4 kcal/mol for PA(H 2 ) (after all, both unpro- 
tonated and protonated systems should be treatable essentially exactly at this level) is in 
part due to the harmonic approximation for the zero-point energy. We have calculated 
CCSD(T)/cc-pVQZ quartic force fields for H 2 and H3 , and found that the anharmonic zero- 
point energies thus obtained lead to PA(H 2 )=101.10 kcal/mol, or only 0.2 kcal/mol higher 
than the evaluated Hunter and Lias value. 

V. CONCLUSIONS 

We have assessed the performance of two recently developed methods for benchmark- 
quality computational thermochemistry, Wl and W2 theory, for a fairly large set of heats 
of formation, as well as for the G2-1 and G2-2 sets of ionization potentials and electron 
affinities, and a number of proton affinities. 

For molecules which exhibit slow basis set convergence, the numerical stability of the Wl 
method is considerably enhanced by substituting the three-point geometric extrapolation of 
the SCF component, E SC f{L) = E S cf,oo + A/B L , by a two-point extrapolation E SCF (L) = 
Escf,oo + A/L 5 which does not involve results with the smallest of the three valence basis 
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sets. 

Wl theory performs excellently for ionization potentials, achieving a mean absolute de- 
viation (MAD) of 0.013 eV for the G2-1 set and 0.018 eV for the G2-2 set (minus CN, B 2 H 4 , 
sec-C3H 7 , Si2H 6 , and CH 3 OF). Both mean and maximum errors are several times smaller 
than those of other, less expensive, theoretical thermochemistry methods. The vastly more 
expensive W2 method (which uses the same correlation methods but larger basis sets) yields 
only a marginal improvement over the Wl method: it appears that the main impediments 
to greater accuracy for these properties are the residual imperfections (higher-order T 3 , T 4 ) 
in the CCSD(T) method. 

The performance of Wl theory for electron affinities is similar, with a MAD of 0.016 
eV for the G2-1 and 0.019 eV for the G2-2 set. The latter value is reduced to 0.016 eV if 
two strongly multireference systems (0 3 and CH 2 NC) are eliminated. For the G2-1 set, W2 
theory (MAD=0.012 eV) does represent a minor improvement over Wl theory, reflecting the 
stronger basis set sensitivity of electron affinities. 

Inner-shell correlation is somewhat more important for ionization potentials than for 
electron affinities: scalar relativistic effects cannot be neglected for either property if results 
of the highest accuracy are desired. For IPs or EAs involving systems in degenerate states, 
spin-orbit splitting corrections are essential. With the exception of CIO (where a fairly large 
active space is required for good results), sufficiently accurate spin-orbit splittings can be 
computed at the CISD/MTsmall level provided the (2s2p) orbitals on second-row atoms are 
included in the correlation. 

Comparison with experiment for the heats of formation in the G2-1 and (part of) the 
G2-2 set is complicated somewhat by the uncertainties in the experimental values: the 
MADs of both Wl and W2 theory are lower than the average experimental uncertainty. For 
parametrizing more approximate methods, Wl and especially W2 level heats of formation 
may well be more suitable than the experimental data. 

Computed proton affinities at the Wl level appear to be converged with the level of 
theory, and agree excellently with experiment. 
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As a final conclusion, we have established that the previously proposed Wl and W2 
theories are in fact valuable and powerful tools for accurate ab initio thermochemistry, with 
mean and maximum absolute errors that are several times smaller than those of more popular 
(and less expensive) schemes such as G2, G3, and CBS-Q/CBS-QB3. 
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SUPPLEMENTARY MATERIAL 

Calculated total energies and geometries of the species discussed in the paper, as well as 
aug-cc-pVnZ basis sets for Li, Be, Na, and Mg, are available on the World Wide Web at the 
Uniform Resource Locator (URL) |http : //theochem. weizmann. ac . il/web/papers/wlw2 
html] as well as in Ref. [64[ . 
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TABLES 



TABLE I. Comparison of different extrapolation procedures for the SCF component (kcal/mol) 



Species 


WI 




Wl' 




W2 






{D,T,Q} 


{T,Q} 


{D,T,Q} 


{T,Q} 


{T,Q,5} 


{Q,5} 




Total Atomization Enerf 


;y 








LiF 


90.05 


93.76 


90.05 


93.76 


93.53 


93.51 


BeH 


50.76 


50.44 


50.76 


50.44 


50.32 


50.34 


SiH 2 ( 3 Bi) 


108.49 


108.29 


109.10 


108.33 


108.25 


108.26 


SiH 3 


182.86 


182.66 


183.72 


182.73 


182.55 


182.58 


SiH 4 


260.24 


259.95 


261.69 


260.03 


259.83 


259.86 


Si 2 H 6 


424.69 


424.37 


426.24 


424.50 


424.14 


424.19 


so 2 


121.51 


121.74 


121.87 


122.00 


121.94 


121.98 


BF 3 


374.54 


374.29 


374.54 


374.29 


374.59 


374.58 


CeHg 


1045.53 


1045.15 


1045.53 


1045.15 






CH3COCH3 


736.37 


736.12 


736.37 


736.12 






CH30CH3 


597.05 


596.85 


597.05 


596.85 






C2H40(oxirane) 


470.00 


469.81 


470.00 


469.81 






S0 3 (a) 


159.50 


159.85 


159.93 


160.20 


159.90 


159.97 


Ionization potential 


Al 


126.73 


126.91 


126.72 


126.91 


126.87 


126.91 


S1H4 


235.45 


235.19 


236.70 


235.22 


235.32 


235.32 


H 2 S( 2 A!) 


270.91 


270.89 


271.02 


270.89 


270.83 


270.81 


C1F 


274.18 


274.08 


274.17 


274.06 


274.15 


274.24 


He 


540.86 


540.62 


540.86 


540.62 


540.71 


540.62 


CF 2 


236.27 


236.16 


236.27 


236.16 






SiH 3 


173.79 


173.72 


174.02 


173.72 






Si 2 H 6 


198.47 


198.17 


199.53 


198.18 






Electron affinity 


Si 


21.92 


22.03 


21.78 


22.03 


22.04 


22.03 


SiH 3 


3.96 


4.04 


3.74 


4.03 


4.31 


4.31 


PH 


-2.54 


-2.53 


-2.59 


-2.53 


-2.50 


-2.50 


HS 


28.06 


28.14 


27.92 


28.13 


28.12 


28.12 


Gl 2 


43.28 


43.16 


43.27 


43.14 


43.21 


43.22 


so 2 


13.93 


14.14 


13.88 


13.79 






(a) Not part of G2-1 


or G2-2 set. 


Raw data 


taken from Ref. [[J. 
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TABLE II. Components of Wl computed electron affinities (kcal/mol) of the G2-1 ion test set 
of molecules. 









( L ) 


Gore 


spm-orbit 


Scalar rel. 


r mai 


7PV17 
Zjr V izj 


vv 1 


Species 


limit 


limit 


limit 


corr. 


splitting 


effects 


Energy 




ri/A 


r\ 

L- 


12.66 


14.36 


1.73 


0.28 


-0.08 


-0.07 


28.88 




28.88 


U 


-12.41 


42.03 


3.81 


0.14 


-0.05 


-0.17 


33.35 




33.35 


T? 

r 


30.19 


44.76 


4.17 


0.17 


-0.38 


-0.26 


78.64 




78.64 


hi 


22.03 


9.10 


1.52 


-0.24 


-0.42 


-0.21 


31.78 




31.78 


r> 

Jr 


-10.59 


25.69 


1.85 


-0.06 


0.28 


-0.24 


16.93 




16.93 




20.84 


25.34 


2.15 


0.00 


-0.09 


-0.29 


47.95 




47.95 


CI 


58.34 


24.18 


2.30 


0.03 


-0.84 


-0.34 


83.66 




83.66 


/—ITT 
OH 


8.71 


16.50 


1.95 


0.25 


-0.04 


-0.06 


27.32 


-0.49 


27.81 


L-J12 


-22.70 


33.34 


3.04 


-0.20 


0.00 


0.02 


13.49 


-1.29 


14.78 


/-ITT 

0H3 


-32.71 


29.24 


3.22 


-0.16 


0.00 


0.00 


-0.40 


-1.06 


0.66 


MTU 


-32.56 


36.54 


3.60 


0.09 


0.08 


-0.09 


7.66 


-0.36 


8.03 


MTU 


OO 1 A 


34.94 


3.99 


0.07 


0.00 


-0.07 


16.81 


-0.63 


1 1 .44 


ATT 

Url 


-3.05 


40.75 


4.40 


0.11 


-0.19 


-0. 14 


41.90 


-0.16 


42.06 


g;tu 
blrl 


16.81 


10.62 


1.51 


-0.17 


-0.20 


-0.19 


28.38 


-0.30 


28.68 


SiH 2 


10.74 


12.21 


1.56 


-0.09 


0.00 


-0.18 


24.25 


-0.66 


24.90 


SiH 3 


4.04 


24.64 


1.99 


0.27 


0.00 


0.11 


31.04 


-0.94 


31.98 


PH 


-2.53 


23.71 


2.13 


-0.05 


0.29 


-0.22 


23.31 


-0.21 


23.52 


PH 2 


4.50 


22.08 


2.29 


-0.04 


0.00 


-0.21 


28.61 


-0.44 


29.05 


HS 


28.14 


23.61 


2.47 


-0.01 


-0.54 


-0.27 


53.41 


-0.09 


53.50 


o 2 


-15.48 


22.84 


2.07 


0.02 


0.23 


-0.15 


9.53 


-0.64 


10.16 


NO 


-14.46 


13.99 


0.68 


0.07 


-0.17 


-0.17 


-0.05 


-0.76 


0.72 


CN 


76.63 


13.10 


-0.16 


0.21 


0.00 


-0.03 


89.75 


-0.01 


89.76 


PO 


14.98 


10.23 


0.46 


-0.15 


-0.33 


-0.30 


24.89 


-0.32 


25.20 


s 2 


22.66 


15.40 


0.57 


-0.16 


0.57 


-0.32 


38.72 


-0.21 


38.93 


Cl 2 


43.16 


11.39 


0.71 


-0.09 


0.00 


-0.16 


55.00 


-0.48 


55.48 
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TABLE III. Deviation of electron affinities (eV) from experiment for the G2-1 test set. 



Expt. a 



Deviation(experiment — theory) 



Species 





EA 


±(uncert.) 


Wl 


W2 


G2 b 


G3 6 


CBS-QB3 C 


c 


1.2629 d 


0.0003 


0.011 


0.007 


0.070 


0.070 


0.082 


o 


1.461122 d 


0.000003 


0.015 


0.012 


0.060 


0.126 


0.087 


F 


3.401190 d 


0.000004 


-0.009 


-0.002 


-0.080 


0.002 


0.035 


Si 


1.38946 e 


0.00006 


0.011 


0.010 


0.030 


0.011 


0.039 


P 


0.7465 d 


0.0003 


0.012 


0.015 


0.110 


0.035 


0.030 


S 


2.077104 d 


0.000001 


-0.002 


0.008 


0.070 


0.013 


-0.017 


CI 


3.61269 d 


0.00006 


-0.015 


0.002 


0.010 


0.007 


-0.065 


CH 


1.238 


0.0078 


0.032 


0.029 


0.110 


0.059 


0.108 


CH 2 


0.652 


0.006 


0.011 


0.002 


-0.010 


0.071 


0.030 


CH 3 


0.08 


0.03 


0.051 


0.034 


0.040 


0.119 


0.091 


NH 


0.37 


0.004 


0.022 


0.008 


0.100 


0.175 


0.108 


NH 2 


0.776 


0.037 


0.020 


0.007 


0.000 


0.078 


0.056 


OH 


1.8277 


0.000044 


0.004 


-0.001 


-0.040 


0.050 


0.061 


SiH 


1.2771 


0.0087 


0.034 


0.031 


0.090 


0.007 


0.082 


SiH2 


1.123 


0.022 


0.043 


0.039 


0.140 


0.048 


0.108 


SiH 3 


1.406 


0.014 


0.019 


0.011 


-0.010 


-0.021 


0.043 


PH 


1.028 


0.01 


0.008 


0.010 


0.070 


0.048 


0.026 


PH 2 


1.271 


0.01 


0.011 


0.013 


0.020 


0.000 


0.004 


HS 


2.317 


0.002 


-0.003 


0.008 


0.060 


-0.003 


-0.013 


2 


0.451 


0.007 


0.010 


-0.003 


-0.030 


0.052 


-0.009 


NO 


0.026 


0.005 


-0.005 


-0.001 


0.090 


0.030 


0.017 


CN 


3.862 


0.005 


-0.031 


-0.026 


-0.110 


-0.067 


-0.048 


PO 


1.092 


0.01 


-0.001 


-0.002 


0.050 


-0.057 


0.039 


s 2 


1.663 


0.04 


-0.025 


-0.018 


0.010 


-0.006 


-0.026 


Cl 2 


2.4 


0.2 


-0.006 


0.004 


0.010 


-0.067 


-0.121 


Mean abs. Err 




0.017 


0.016 


0.012 


0.057 


0.049 


0.054 


Max. abs. Err 




0.200 


0.051 


0.039 


0.140 


0.175 


0.121 



a Unless otherwise indicated, experimental values are those from [35] 
b G2 and G3 values from §] 
c CBS-QB3 values from § 

d CRC Handbook of Chemistry and Physics, 78th ed. (CRC, Boca Raton, FL, 1997). 

6 J. Thogersen, L.D. Steele, M. Scheer, C.A. Brodie, and H.K. Haugen, J. Phys. B 29, 1323 

(1996) 
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TABLE IV. Calculated and experimental spin-orbit contributions (eV). 



Calc. spin-orbit splitting 



Spin-orbit splitting 
from experimental 





HF 


CISD a 


CISD+subval 6 


experimental fine 


OH 


0.00848 


0.00833 


0.00838 


0.00863 


CH 


0.00165 


0.00162 


0.00164 


0.00173 


SH 


0.02177 


0.02104 


0.02360 


0.02337 


NO 


0.00692 


0.00729 


0.00732 


0.00743 


OF 


0.01067 


0.01117 


0.01117 


— 


C10 c 


0.01822 


0.01801 


0.01917 


0.01971 


CIO 


0.01223 


0.01391 


0.01378 




NCCN+ 


0.00305 


0.00333 


0.00333 


— 


CS+ 


0.02343 


0.02437 


0.02678 


— 


OCS+ 


0.02198 


0.02157 


0.02381 


— 


CO+ 


0.00876 


0.00936 


0.00939 


— 


N+( 2 n) 


0.00415 


0.00439 


0.00439 


0.00463 


HCCH+ 


0.00172 


0.00183 


0.00183 


— 


NH+ 


0.00480 


0.00475 


0.00480 


0.00482 


PH+ 


0.01678 


0.01605 


0.01852 


— 


C1F+ 


0.03911 


0.03810 


0.04174 


0.03906 


HF+ 


0.01761 


0.01732 


0.01740 


0.01815 


HC1+ 


0.03742 


0.03627 


0.03999 


0.04018 


C1+ 


0.03965 


0.04031 


0.04416 


0.03998 


ot 


0.01181 


0.01199 


0.01199 


0.01223 


pi 


0.01328 


0.01365 


0.01521 


0.01612 


o + 
D 2 


0.02649 


0.02649 


0.02970 


0.02914 


NH" 


0.00344 


0.00341 


0.00344 




o 2 - 


0.00941 


0.00976 


0.00976 


0.00992 


SiH 


0.00787 


0.00748 


0.00882 


0.00885 


PO 


0.01293 


0.01275 


0.01426 


0.01389 


PH~ 


0.01117 


0.01086 


0.01241 






0.02169 


0.02214 


0.02455 


(0.026) 



62 



Experimental fine structure data from Ref. 
a CISD with valence correlation only 

b CISD with valence correlation, and 2s2p correlation in second-row atoms 
c CASSCF calculated values 



32 



TABLE V. Deviation of electron affinities (eV) from experiment for the G2-2 test set. 



Expt. a 



Deviation(cxpcrimcnt — theory) 



Species 
























EA 


±(uncert.) 


Wl 


G2 b 


G3 6 


Li 


0.61759 


0.00022 


0.001 c 


-0.132 


-0.124 


B 


0.27972 d 


0.00003 


0.021 


0.090 


0.076 


Na 


0.547951 


0.000044 


-0.002 c 


-0.132 


-0.159 


Al 


0.43283 6 


0.00005 


0.020 


0.083 


0.043 


C 2 


3.273 


0.008 


0.010 


0.173 


0.116 


C2O 


2.289 


0.018 


-0.016 


-0.041 


-0.001 


CF 2 


0.179 


0.005 


0.019 


0.089 


0.001 


NCO 


3.609 


0.005 


-0.021 


-0.011 


0.032 


NO2 


2.273 


0.005 


-0.014 


-0.067 


-0.008 


O3 


2.103 


0.004 


-0.066 


0.033 


0.000 


OF 


2.272 


0.006 


-0.041 


-0.028 


0.021 


SO2 


1.107 


0.008 


-0.029 


-0.053 


-0.077 


S 2 


1.877 


0.008 


-0.026 


-0.043 


-0.113 


CCH 


2.969 


0.006 


0.007 


-0.121 


-0.027 


C2H3 


0.667 


0.024 


-0.007 


-0.083 


0.012 


CH 2 CC 


1.794 


0.008 


-0.024 


0.054 


-0.027 


CH 2 CCH 


0.893 


0.005 


-0.016 


-0.097 


-0.013 


CH2CHCH2 


0.481 


0.008 


-0.005 


-0.039 


0.039 


HCO 


0.313 


0.005 


0.000 


-0.027 


0.005 


HCF 


0.542 


0.005 


0.004 


0.082 


0.013 


CH3O 


1.570 


0.005 


0.014 


-0.050 


0.017 


CH 3 S 


1.861 


0.004 


-0.010 


-0.009 


0.001 


CH 2 S 


0.465 


0.023 


-0.041 


0.075 


0.001 


CH 2 CN 


1.544/ 


0.006 


-0.002 


-0.036 


0.026 


CH 2 NC 


1.058-f 


0.026 


-0.075 


-0.122 


-0.048 


CHCO 


2.35 


0.022 


0.008 


-0.010 


0.039 


CH2CHO 


1.824157 


0.000044 


-0.008 


-0.046 


-0.010 


CH3CO 


0.423 


0.037 


0.027 


-0.017 


0.020 


CH3CH2O 


1.713 


0.005 


-0.047 


-0.097 


-0.039 


LiH 


0.342 


0.012 


0.013 


0.022 


-0.053 


HNO 


0.338 


0.015 


0.003 


0.088 


0.043 


H0 2 


1.078 


0.017 


-0.002 


-0.032 


0.029 


Mean abs. Err 




0.010 


0.019 


0.065 


0.039 


Max. abs. Err 




0.037 


0.075 


0.173 


0.159 



a Unless otherwise indicated, experimental values are those from |55] 

b G2 and G3 values from |fl 

c augmented basis sets were used 

d M. Scheer, R.C. Bilodeau, and H.K. Haugen, Phys. Rev. Lett. 80, 2562 (1998) 

e M. Scheer, R.C. Bilodeau, J. Thogersen, and H.K. Haugen, Phys. Rev. A 57, 1493 (1998) 

f Ref. 1721 
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TABLE VI. Components of Wl computed ionization potentials (kcal/mol) of the G2-1 ion test 



set of molecules. 





SCF 


CCSD 


(T) 


Core 


spin-orbit 


Scalar rel. 


Final 


ZPVE 


Wl 


Species 


limit 


limit 


limit 


corr. 


splitting 


effects 


Fnorp'v 




IP 


Li 


123.18 


0.00 


0.00 


1.02 


0.00 


0.01 


124.21 




124.21 


Be 


185.51 


29.03 


0.00 


0.40 


0.00 


0.02 


214.97 




214.97 


13 


182.94 


6.08 


1.05 


0.88 


0.03 


-0.05 


190.93 




190.93 




248.73 


9.29 


0.74 


0.76 


-0.04 


-0.09 


259.40 




259.40 




321.79 


12.29 


0.67 


0.67 


-0.23 


-0.14 


OOO . uu 




OOO .uu 


o 


275.90 


36.06 


1.30 


0.33 


0.22 


-0.20 


313.62 




313.62 


F 


362.02 


37.92 


1.61 


0.34 


-0.09 


-0.29 


401.51 




401.51 


Na 


114.18 


0.00 


0.00 


3.95 


0.00 


0.15 


118.28 




118.28 




152.37 


21.36 


0.00 


2.07 


0.00 


0.22 


176.02 




176.02 


Al 


126.91 


9.93 


0.97 


-0.33 


0.21 


-0.22 


137.48 




137.48 


Si 


176.45 


10.32 


1.14 


-0.08 


-0.09 


-0.25 


187.49 




187.49 


p 


231.39 


9.90 


1.24 


0.08 


-0.85 


-0.28 


241.48 




241.48 


s 


209 63 


27.22 


1.22 


0.23 


0.56 


-0.32 


238.54 




238.54 


CI 


271.73 


26.08 


1.36 


0.25 


-0.13 


-0.36 


298.94 




298.94 


CH^ (T)nj ration 1 


278.37 


18.94 


0.68 


0.14 


0.00 


-0.03 


298.10 


5.43 


292.67 


CHa(C,o ration ) a 


271.66 


22.52 


1.05 


0.27 


0.00 


-0.03 


295.47 


3.93 


291.54 




202.33 


31.12 


2.35 


0.09 


0.00 


-0.02 


235.87 


1.04 


234.83 


OH 


264.15 


34.73 


1.88 


0.28 


0.19 


-0.16 


301.07 


0.93 


300.14 


H 2 


255 86 


34.37 


2.39 


0.24 


0.00 


-0.13 


292.73 


1.68 


291.05 


FH 


331.78 


38.03 


2.28 


0.27 


0.40 


-0.24 


371.72 


1.54 


370.18 


g;H 4 


235.19 


19.71 


0.61 


0.23 


0.00 


-0.25 


255.49 


1.71 


253.78 


PH 


221.61 


11.75 


1.25 


0.17 


-0.43 


-0.25 


234.11 


-0.01 


234.12 


PH 2 


211.53 


13.52 


1.33 


0.26 


0.00 


-0.23 


226.42 


-0.03 


226.45 




198.75 


26.57 


1.89 


0.49 


0.00 


0.19 


227.89 


0.21 


227.67 


HS 


213.11 


25.26 


1.54 


0.21 


0.54 


-0.29 


240.37 


0.20 


240.17 


HoSC 2 Bi "l 


216.34 


23.62 


1.80 


0.19 


0.00 


-0.27 


241.67 


0.39 


241.28 




270.89 


22.49 


1.40 


0.30 


0.00 


0.24 


295.31 


0.39 


294.91 


HC1 


269.43 


24.27 


1.71 


0.21 


-0.92 


-0.33 


294.37 


0.45 


293.92 




227.74 


31.98 


3.30 


0.64 


-0.04 


-0.04 


263 59 


0.61 


262.98 


C2H4 


207.05 


33.10 


3.40 


0.49 


0.00 


-0.04 


244.01 


1.44 


242.56 


CO 


306.44 


16.79 


-0.20 


0.22 


0.00 


0.02 


323.27 


-0.11 


323.38 


N 2 ( 2 S+ cation) 
N 2 ( 2 ncation) 


368.35 


-4.69 


-3.62 


0.56 


0.00 


-0.07 


360.52 


0.17 


360.36 


354.09 


29.27 


1.98 


0.96 


-0.10 


-0.16 


386.03 


0.17 


385.87 


O2 


271.54 


7.58 


-0.81 


0.32 


-0.28 


-0.21 


278.15 


-0.59 


278.75 


P 2 


214.65 


24.92 


3.35 


0.21 


-0.35 


-0.36 


242.42 


0.06 


242.36 


S2 


208.84 


8.59 


-0.61 


0.04 


-0.68 


-0.39 


215.77 


-0.17 


215.94 


Cl 2 


255.59 


11.60 


-0.79 


0.09 


-1.02 


-0.40 


265.07 


-0.13 


265.20 


C1F 


274.08 


18.66 


0.13 


0.19 


-0.96 


-0.32 


291.78 


-0.20 


291.98 


CS 


243.16 


17.48 


0.91 


0.01 


0.00 


-0.06 


261.51 


-0.15 


261.66 



a geometry optimized at CCSD(T)/cc-pVTZ level. B3LYP/cc-pVTZ optimization erro- 
neously yields D 2 d structure (see text and Table VIII). 
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TABLE VII. Deviation of ionization potentials (eV) from experiment for the G2-1 test set. 



Expt. a Deviation(cxporimcnt — theory) 
Species 





IP 


±(uncert.) 


Wl 


W2 


G2 b 


G3 6 


CBS-QB3 C 


Li 


5.39172 


0.00001 


0.005 


0.005 


0.082 


-0.007 


0.084 


Be 


9.32263 


0.00001 


0.001 


0.005 


-0.090 


-0.135 


-0.048 


13 


8.29802 


0.00002 


0.019 


0.007 


0.100 


0.063 


0.074 


C 


11.2603 


0.0001 


0.012 


0.010 


0.080 


0.051 


0.070 


N 


14.534 


0.001 


0.004 


0.000 


0.060 


0.029 


0.056 





13.618 


0.001 


0.018 


0.005 


0.080 


0.071 


0.022 


F 


17.423 


0.001 


0.012 


0.002 


0.030 


0.034 


-0.052 


Na 


5.139 


0.001 


0.010 


0.010 


0.189 


0.027 


0.009 


Mg 


7.646 


0.001 


0.013 


0.013 


-0.004 


-0.137 


0.053 


Al 


5.986 


0.001 


0.024 


0.023 


0.050 


0.028 


0.061 


Si 


8.15166 


0.00003 


0.021 


0.018 


0.050 


0.025 


0.060 


P 


10.48669 


0.00001 


0.015 


0.011 


0.037 


0.023 


0.079 


S 


10.360 


0.001 


0.016 


0.014 


0.160 


0.092 


0.091 


CI 


12.968 


0.001 


0.005 


0.007 


0.120 


0.076 


0.047 


CH4(D2d cation) 


12.61 


0.01 


-0.081 


-0.082 


-0.060 


-0.043 


-0.100 


CH 4 (C2i, cation) d 


12.61 


0.01 


-0.032 


-0.033 


-0.060 


-0.043 


-0.100 


NH 3 


10.18 


0.09 


-0.003 


-0.004 


0.101 


0.042 


0.115 


OH 


13.017 


0.0002 


0.002 


0.001 


0.030 


0.082 


-0.013 


H 2 


12.6223 


0.0003 


0.001 


0.006 


-0.010 


0.030 


-0.039 


FH 


16.044 


0.003 


-0.008 


-0.016 


-0.050 


0.000 


-0.091 


SiH 4 


11.0 


0.02 


-0.005 


0.006 


-0.010 


-0.023 


-0.009 


PH 


10.149 


0.008 


-0.003 


-0.006 


0.060 


-0.037 


0.043 


PH2 


9.824 


0.002 


0.004 


0.003 


0.100 


0.007 


0.061 


PH 3 


9.87 


0.002 


-0.003 


-0.006 


0.000 


-0.012 


0.026 


HS 


10.4219 


0.0004 


0.007 


0.007 


0.060 


0.097 


0.026 


H 2 S( 2 Bi) 


10.453 


0.001 


-0.010 


-0.008 


0.040 


0.011 


0.022 


H 2 S( 2 Ai) 


12.76 


0.036 


-0.029 


-0.029 


0.030 


-0.002 


-0.004 


HC1 


12.747 


0.002 


0.001 


0.008 


0.030 


0.029 


-0.026 


C2H2 


11.403 


0.0003 


-0.001 


-0.004 


-0.020 


-0.006 


-0.039 


C2H4 


10.5138 


0.0006 


-0.005 


-0.001 


-0.070 


-0.045 


-0.035 


CO 


14.0142 


0.0003 


-0.009 


-0.014 


0.000 


-0.001 


-0.056 


N 2 ( 2 £+cation) 


15.581 


0.008 


-0.046 


-0.046 


0.020 


0.018 


-0.030 


N 2 ( 2 ncation) 


16.699 


0.001 


-0.034 


-0.049 


0.030 


0.030 


-0.052 


O2 


12.0697 


0.0002 


-0.018 


-0.024 


-0.100 


-0.176 


-0.095 


P 2 


10.567 


0.002 


0.057 


0.047 


-0.010 


0.017 


-0.013 


S 2 


9.356 


0.002 


-0.008 


-0.011 


0.080 


-0.023 


-0.048 


Cl 2 


11.481 


0.003 


-0.019 


-0.008 


-0.010 


-0.045 


-0.039 


C1F 


12.66 


0.01 


-0.002 


0.005 


-0.070 


-0.002 


-0.110 


CS 


11.33 


0.01 


-0.017 


-0.017 


-0.090 


-0.061 


0.048 


Mean abs. error. 
Max abs. error. 




0.006 
0.090 


0.013 
0.057 


0.013 
0.049 


0.058 
0.189 


0.043 
0.176 


0.051 
0.115 



a Experimental values from |85] 
b G2 and G3 values from £| 
c CBS-QB3 values from § 

d CCSD(T)/cc-pVTZ geometry. B3LYP/cc-pVTZ optimization erroneously yields Did struc- 
ture (see text and Table VIII). 
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TABLE VIII. Equilibrium structure (A, degrees) and harmonic frequencies (cm 1 ) of the Ci v 
structure of CH4" obtained with the cc-pVTZ basis set and a variety of electronic structure methods 





CCSD(T) 


B3LYP 


BLYP 


BHLYP 


BHPW91 


B3PW91 


B3P86 


mPWlPW91 


mPWIK 


rem 


1.0826 


1.0827 


1.0888 


1.0752 


1.0764 


1.0838 


1.0829 


1.0825 


1.0781 


rcH2 


1.1864 


1.1867 


1.1955 


1.1766 


1.1769 


1.1869 


1.1858 


1.1851 


1.1790 


6hichi> 


125.62 


124.20 


124.12 


124.23 


124.88 


124.75 


124.72 


124.80 


124.84 




55.01 


57.56 


58.22 


57.00 


55.76 


56.52 


56.58 


56.34 


55.97 


B 2 


3308 


3281 


3208 


3370 


3384 


3292 


3298 


3308 


3362 


A! 


3165 


3147 


3078 


3232 


3242 


3154 


3161 


3170 


3221 


A 1 


2572 


2534 


2460 


2620 


2641 


2552 


2557 


2569 


2621 


B 2 


2292 


2188 


2097 


2281 


2353 


2246 


2248 


2267 


2328 


A x 


1585 


1545 


1507 


1593 


1615 


1562 


1564 


1572 


1602 


Bi 


1311 


1296 


1269 


1337 


1331 


1292 


1292 


1297 


1321 


Ai 


1201 


1159 


1126 


1184 


1226 


1187 


1191 


1196 


1218 


B 2 


894 


901 


875 


933 


912 


885 


889 


888 


906 


A 2 


469 


227i 


417? 


358 


339 


233i 


256i 


156i 


269 
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TABLE IX. Deviation of ionization potentials (eV) from experiment for the G2-2 test set. 



Expt. a 



Deviation(experiment — theory) 



Species 





IP 


±(uncert.) 


Wl 


G2 6 


G3 b 


H 


13.599 


0.001 


-0.007 


0.004 


0.035 


He 


24.588 


0.001 


0.000 


-0.048 


-0.018 


Ne 


21.565 


0.001 


-0.012 


0.048 


0.013 


Ar 


15.760 


0.001 


-0.005 


-0.067 


-0.071 


BF 3 


15.71 


0.1 


-0.028 


-0.099 


-0.103 


CO2 


13.778 


0.002 


-0.039 


-0.071 


-0.079 


CF2 


11.42 


0.01 


-0.010 


0.024 


0.019 


OCS 


11.185 


0.002 


0.006 


-0.023 


-0.006 


CS 2 


10.08 


0.002 


-0.003 


0.020 


0.015 


CH 2 


10.396 


0.003 


0.022 


-0.084 


-0.002 


CH 3 


9.837 


0.005 


-0.003 


-0.054 


0.033 


C 2 H 5 ( 2 A') 


8.117 


0.008 


-0.009 


-0.047 


0.048 


C3H4 (cyclopropene) 


9.668 


0.005 


0.032 


0.093 


0.067 


CH 2 CCH 2 


9.6878 


0.002 


-0.014 


0.047 


0.008 


sec-C3H7 


7.37 


0.02 


-0.079 


0.028 


0.093 


C 6 H 6 


9.24384 


0.00006 


-0.001 


0.084 


0.066 


CN( 1 S+ cation) c 


14.03 


0.02 


0.165 


0.297 


0.171 


CN( 3 II cation) c 


14.03 


0.02 


-0.039 






CHO 


8.14 


0.04 


-0.010 


■ 




CHaOH^A) 


7.56 


0.01 


0.016 


-0.110 


- 


CH 3 0( 2 A') 


10.72 


0.01 


-0.006 


0.060 


0.021 


CH3OH 


10.85 


0.03 


-0.035 


0.116 


0.077 


CH 3 F 


12.54 


0.01 


-0.021 


0.161 


0.144 


CH 2 S 


9.376 


0.003 


-0.013 


-0.001 


-0.018 


CH 2 SH 


7.536 


0.003 


-0.026 


-0.116 


-0.034 


CH3SH 


9.443 


0.002 


-0.007 


0.019 


0.015 


CH3CI 


11.265 


0.003 


-0.022 


0.036 


0.027 


CH3CHO 


10.227 


0.005 


-0.031 


0.081 


0.046 


CH3OF 


11.34 


0.008 


-0.078 


0.065 


0.060 


C 2 R*4S(thiirane) 


9.051 


0.006 


-0.013 


0.016 


0.016 


NCCN 


13.374 


0.008 


-0.046 


0.017 


0.026 


C4H 4 0(furan) 


8.91 


0.01 


0.032 


-0.003 


-0.020 


B 2 H4 


9.70 


0.02 


0.129 


0.128 


-0.091 


NH 


13.49 


0.01 


0.019 


-0.077 


0.005 


NH 2 


11.14 


0.01 


-0.033 


0.035 


-0.021 


N 2 H 2 


9.589 


0.007 


-0.010 


0.129 


0.086 


N 2 H 3 


7.61 


0.01 


-0.005 


-0.069 


0.005 


HOF 


12.71 


0.01 


-0.022 


-0.004 


0.009 


SiHaCAO 


9.15 


0.02 


0.015 


0.022 


0.043 


SiH 3 


8.135 


0.005 


0.022 


-0.082 


0.017 


Si 2 H 2 


8.2 


0.02 


0.008 


0.078 


0.057 


Si 2 H 4 


8.09 


0.03 


0.030 


0.032 


0.054 


Si 2 H 6 


9.74 


0.02 


0.091 e 


-0.035 


-0.044 


Mean abs. Err 




0.011 


0.025 


0.062 


0.044 


w/o B 2 H4 




0.011 


0.022 


0.062 


0.043 


w/o B 2 H 4 , sec-C 3 H 7 , CH3OF, Si 2 H 6 




0.010 


0.017 


0.063 


0.041 


Max abs. Err 




0.100 


0.129 


0.161 


0.144 


w/o B 2 H 4 , sec-C 3 H 7 , CH3OF, Si 2 H 6 




0.100 


0.046 


0.161 


0.144 



a Experimental values from [85] 
b G2 and G3 values from J§ 

c Ionization potential of 13.60 eV given in Ref. |^{| is propagated transcription error (see text). Measured 
value of 14.03±0.02 eV ]73] corresponds to unspecified state of cation: may be 3 II for symmetry reasons. 
Consequently, CN has been excluded from the error statistics for all methods. 
d W2 theory: IP=9.57 eV, experiment-theory=0.127 eV 
e W2 theory: IP=9.66 eV, experiment-theory=0.084 eV 
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TABLE X. Deviation of heats of formation (kcal/mol) from experiment for the G2-1 test set. 



Expt. a Deviation(experimcnt — theory) 
Species 







zb(uncert.) 


Ref. 


Wl 


W2 


G2 6 


G3 6 


CBS-Q C 


T ill 

J_j1x1 


00. Dl 


n m 
U.Ul 


T 

J 


n a 
u. 1 


n q 
U . 


n 
u.y 


n fi 

U.O 


n 1 
-U.l 
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ol. f U 


1 nn 
l.UU 


1_[ 
11 
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U.o 


1 c; 
-1.0 


n c; 
-U.O 
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-U.o 
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VT 
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u.i 
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U.l 
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U.Z 
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VI 
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oU. ( 4 


n H7 
U.U / 


c 

V. 


n q 
U.o 


n ("i 
u.y 


Q O 

-o.Z 


n 

-u.y 


"i 7, 
l.o 


v.'i2 


0.00 


0.00 


s 


-0.6 


0.2 


-1.4 


-1.1 


1.7 


NaCl 




U.UU 


T 

J 




-u.u 


1 2 


1 4 


\>.o 


DlO 


o/i nn 

-Z4.UU 


m 

Z.Ul 


J 


n 
-u.y 


n fi 
-U.O 


-1.1 


n 1 

-U. 1 


"1 fi 
1.0 


cs 


66.86 


0.23 


H 


-0.9 


-0.5 


1.0 


1.1 


1.3 


so 


1.20 


0.31 


J 


-0.4 


0.1 


-2.6 


-0.5 


0.7 


CIO 


24.19 


0.50 


J 


-1.3 


-0.5 


-2.2 


-1.7 


-0.3 


FC1 


-13.20 


0.10 


J 


0.0 


0.0 


0.7 


-0.7 


0.5 


Si 2 H 6 


19.05 


0.31 


L 


2.2 


1.9 


2.8 


1.3 


3.3 


CH3CI 


-19.57 


0.14 


G 


0.7 


0.6 


0.9 


-0.1 


1.1 


CH3SH 


-5.46 


0.14 


GO 


0.8 




-0.2 


-0.4 


0.5 


HOC1 


-17.81 


0.50 


J 


0.5 


0.5 


0.5 


-0.4 


1.1 


SO2 


-70.94 


0.05 


C 


-0.7 


0.2 


-4.9 


-3.8 


-0.3 


Mean abs. error 




0.54 




0.6 


0.5 


1.2 


1.1 


1.0 


Max abs. error 




3.11 




2.2 


1.9 


4.9 


4.0 


3.3 



a Experimental values from: 
J = JANAF H Tables 
H = Huber and Herzberg 
C 
G 



CODATA Values 
Gurvich et al. |8^] compilation 
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S=standard state 

DO = Doncaster, and Walsh, Int. J. Chem. Kinet. 13, 503 (1981) 
L = Lias et al. |88[ compilation 

GO = W.D. Good, J.L. Lacina, J.P. McCullough, J. Phys. Chem. 65, 2229 (1961) 
E = Based on singlet-triplet splitting determined by A.R.W. McKelllar, P.R. Bunker, T.J. 
Sears, K.M. Evenson, R.J. Saykally and S.R. Langoff, J. Chem. Phys. 79, 5251 (1983) 
B = J. Berkowitz, J.P. Greene, H. Cho, and B. Ruscic, J. Chem. Phys. 86, 1235 (1987) 
b G2 and G3 values from 
c CBS-Q values from 
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TABLE XI. Deviation of heats of formation (kcal/mol) from experiment for selected molecules 
from the G2-2 test set. 

Expt. a Deviation(cxporimcnt — theory) 
Species 





AH° 


-1— ( 11 nr , £>i*1~ 1 
1 1 LI 1 1 l x u . t 


Ref. 


Wl 


W2 


G2 & 


G3 b 




r>r 3 


071 X. 

-4 ( 1.0 


U.Z 




l.O 


1 

-U.l 


-U.l 


n ft 
-U.o 


1 
-l.Z 


Dpi 


Oft Q 

-yo.o 


U.O 


T 
■ J 


U.l 


1 

l.U 




z.u 


U.U 


l.U 


<^r 4 


-ZZo.U 


U.Z 


r 1 

V T 


"1 7 




u.z 


. 



U . :f 


. 




-00. y 


U.O 


1^ 
U 


O 1 

u. 1 



U.Z 


I 

i.y 


Z.U 


0. 1 


nc 


z / .y 


U.Z 




n 1 
U. 1 


U. ( 



z.U 


o.Z 


0.4 


N20 


19.5 


0.1 


G 


-0.8 


-1.2 


-0.7 


-1.9 


-0.3 




12.6 


0.1 




-1.6 


-1.6 


1.0 


-0.8 


1. r 




Q1 ft 

-ol.O 


U.O 


.) 


n 1 
U. 1 




6. i 


U.U 


Z.o 


r 2^ 


^ 

o.y 


n 1 
u. 1 


r 1 


-U.o 


-l.U 


U . 


-U.o 


O Q 
U.O 


Crl2-r 2 


-lUo. 1 


U.Z 




U.o 




Z. f 


U.o 


1 Q 
l.O 




177 
1 / . f 


n 1 
u. 1 




U.o 




-U.4 


-U.l 


1 n 
-l.U 


liL-UUn 


on c; 

-yu.o 


n 1 
U. 1 




n ft 
u.o 





Z.U 


i 
U.l 


1. 1 


C2H40(oxirane) 


1 Q ft 

-1Z.0 


n 
U.Z 


■ J 


n 
-U.Z 




1 Q 

1.0 


U.U 




U.Z 




-ou. 1 


n 9 
u. z 


T 


fl 7 




9 O 


Q 
O.J 


1 ^ 

l.U 


CH3OCH3 


-44.0 


0.1 


T 


1.1 




2.0 


0.1 


0.5 


G2H5GI 


-26.8 


0.2 


1 


1.0 




0.8 


-0.1 


0.9 


CH3COCH3 


-51.9 


0.2 


P 


1.1 




1.1 


0.1 


0.0 


HS 


34.2 


0.7 


JB 


0.5 


0.6 


-0.3 


0.5 


0.4 


CCH 


135.9 


1.2 


G 


0.1 


-0.1 


-2.8 


-0.4 


-1.3 


C2H3 


71.5 


1.2 


TS 


0.9 




-1.2 


1.0 


-0.3 


CH3CO 


-2.9 


0.7 


TS 


0.2 




-0.1 


-0.4 


-0.4 


CH2OH 


-4.3 


0.6 


G 


0.4 




-0.5 


-0.4 


-0.5 


CH3CH2O 


-3.7 


0.8 


.JB 


-0.8 




-1.4 


-1.2 


-1.9 


CH 3 S 


29.8 


0.4 


.JB 


1.2 




-0.1 


0.8 


0.9 


C2H5 


28.4 


0.5 


TS 


0.4 




-1.5 


-0.3 


-1.7 


N0 2 


8.2 


0.1 


G 


-0.8 


-1.2 


1.0 


0.1 


2.6 


Mean abs. error 
Max abs. error 




0.4 
1.2 




0.7 
1.7 




1.5 
5.5 


0.7 
3.2 


1.5 
5.1 



a Experimental values from: 

J = JANAF g Tables 

C = CODATA g Values 

G = Gurvich et al. |57J compilation 

T = TRC |K| compilation 

TS = Tsang f9]| compilation 

AN = X. An, M. Mansson, J. Chem. Thermodyn. 15, 287 (1983) 

P = Pedley et al. [|92] compilations 

JB= Values reported by Berkowitz et al. [72 

b G2 and G3 values from || 

c CBS-Q values from |§ 
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TABLE XII. Computed and observed proton affinities (kcai/mol) 





Experiment 




Experiment— theory 






Ref. |4| 


Wl 


W2 


G3 


NH 3 


204.0 


-0.1 


0.1 


0.9 


H 2 


165.2 


0.2 


0.3 


1.8 


C2H2 


153.3 


-0.8 


-0.9 


0.5 


SiH 4 


152.9 


-0.3 


-0.2 


0.6 


PH 3 


187.6 


0.4 


0.5 


2.3 


H 2 S 


168.5 


-0.6 


-0.6 


1.5 


HC1 


133.1 


-0.7 


-0.8 


0.5 


H 2 


100.9 


-0.4 


-0.4° 


1.6 


Mean abs. error 




0.44 


0.49 


1.2 


Max. abs. error 




0.8 


0.9 


2.3 



(a) Using CCSD(T)/cc-pVQZ anharmonic zero point energies: -0.2 kcal/mol. 
(ZPVE(H 3 f )=12.56 rather than 12.31 kcal/mol; ZPVE(H 2 )=6.21 compared to 6.22 
kcal/mol.) 
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S-0 



TABLES 

TABLE S-I. Components of Wl computed electron affinities (kcal/mol) of the G2-2 ion test 



set of molecules. 





bCl 




(T) 


Core 


spin-orbit 


1 . . , . .. . . 1 

scalar rel. 


Final 


ZPVE 


Wl 


Species 


limit 


limit 


limit 


corr. 


splitting 


effects 


Energy 




EA 


Li a 


-2.83 


17.07 


0.00 


0.02 


0.00 


0.00 


14.26 




14.26 


B 


-6.18 


10.40 


1.55 


0.24 


-0.01 


-0.04 


5.96 




5.96 


Na a 


-2.38 


14.96 


0.00 


0.11 


0.00 


0.00 


12.69 




12.69 


Al 


0.95 


7.96 


1.33 


-0.46 


-0.09 


-0.17 


9.52 




9.52 


C2 


99.80 


-18.58 


-6.52 


0.44 


0.00 


0.09 


75.23 


-0.02 


75.26 


C2O 


24.54 


26.33 


2.37 


0.34 


0.00 


-0.13 


53.46 


0.31 


53.15 


CF2 


-13.51 


14.91 


0.97 


0.09 


0.00 


-0.14 


2.32 


-1.37 


3.69 


NCO 


57.63 


24.67 


1.86 


0.26 


0.00 


-0.16 


84.26 


0.54 


83.72 


NO2 


36.83 


15.77 


-0.27 


-0.08 


0.00 


-0.08 


52.18 


-0.56 


52.74 


O3 


54.47 


-0.13 


-5.22 


0.17 


0.00 


-0.17 


49.12 


-0.89 


50.01 


Ol 


19.52 


29.80 


3.39 


0.11 


0.26 


-0.20 


52.89 


-0.45 


53.34 


SO2 


13.88 


12.12 


-0.23 


-0.07 


0.00 


-0.25 


25.44 


-0.76 


26.20 


S2O 


37.09 


8.08 


-1.52 


-0.46 


0.00 


0.16 


43.36 


-0.52 


43.88 


CCH 


44.31 


21.94 


2.53 


-0.18 


0.00 


0.04 


68.64 


0.34 


68.31 


C2H3 


-13.58 


25.81 


2.81 


-0.29 


0.00 


0.00 


14.75 


-0.78 


15.53 


CH2CC 


20.44 


19.51 


0.89 


0.25 


0.00 


-0.08 


41.01 


-0.91 


41.92 


("ITT f If (IT 


-3.43 


21.36 


2.92 


-0.22 


0.00 


-0.02 


20.61 


-0.35 


20.96 


/~1TT /~1TT/— ITT 
CH2CHCH2 


-13.06 


19.60 


2.66 


0.05 


0.00 


-0.06 


9.19 


-2.00 


11.20 


HCO 


-17.58 


20.97 


1.78 


-0.29 


0.00 


-0.03 


4.86 


-2.37 


7.22 


HCF 


-6.34 


16.12 


1.35 


0.11 


0.00 


-0.10 


11.13 


-1.29 


12.42 


CH3O 


-4.98 


36.13 


3.88 


-0.08 


0.00 


-0.14 


34.81 


-1.07 


35.88 


CH 3 S 


17.42 


23.73 


2.50 


-0.13 


0.00 


-0.25 


43.26 


0.12 


43.15 


CH 2 S 


-1.96 


11.96 


0.16 


-0.31 


0.00 


0.01 


9.86 


-1.81 


11.68 


CH 2 CN 


10.39 


22.44 


2.46 


0.04 


0.00 


-0.08 


35.25 


-0.39 


35.64 


CH 2 NC 


-4.68 


27.58 


3.01 


0.01 


0.00 


-0.03 


25.90 


-0.24 


26.14 


CHCO 


27.45 


24.20 


2.42 


0.15 


0.00 


-0.09 


54.14 


0.14 


54.00 


CH2CHO 


12.94 


25.69 


3.13 


0.11 


0.00 


-0.09 


41.78 


-0.48 


42.26 


CH3CO 


-19.31 


24.19 


2.67 


-0.36 


0.00 


0.00 


7.19 


-1.95 


9.14 


CH3CH2O 


-1.89 


37.06 


4.22 


-0.09 


0.00 


-0.12 


39.18 


-1.39 


40.58 


LiH 


5.99 


1.04 


0.20 


0.08 


0.00 


0.00 


7.29 


-0.30 


7.59 


HNO 


-11.26 


16.84 


1.03 


0.04 


0.00 


-0.12 


6.52 


-1.21 


7.74 


H0 2 


-10.01 


30.79 


3.44 


0.02 


0.00 


-0.11 


24.13 


-0.77 


24.90 



a augmented basis sets were used 
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TABLE S— II. Components of Wl computed ionization potentials (kcal/mol) of the G2-2 ion 



test set of molecules. 





SCF 


CCSD 


(T) 


Core 


spin-orbit 


Scalar rel. 


Final 


ZPVE 


Wl 


Species 


limit 


limit 


limit 


corr. 


splitting 


effects 


Energy 




IP 



H 


313. 


.75 





.00 





.00 


0. 


00 





00 


0. 


00 


313 


.75 




313.75 


He 


540. 


62 


26 


39 


0. 


00 





00 





00 


0. 


00 


567 


.02 




567.02 


Ne 


456. 


.35 


39 


.81 


1 


.88 


0, 


34 


-0 


.38 


-0. 


41 


497 


.59 




497.59 


Ar 


339. 


.49 


24. 


.12 


1 


.45 


0. 


23 


-1 


.36 


-0. 


38 


363 


.55 




363.55 


BF 3 


362. 


09 


4 


.55 


-2 


36 


0. 


.41 





.00 


-0. 


34 


364 


.34 


1.42 


362.92 


C0 2 


312. 


.49 


8 


.18 


-0. 


.20 


0. 


40 


-0 


.22 


-0. 


21 


320 


.44 


1.81 


318.63 


CF 2 


236. 


.16 


25. 


.70 


0. 


.91 


-0 


.13 


0. 


00 


0. 


01 


262 


.65 


-0.93 


263.58 


OCS 


233. 


54 


22 


56 


2 


.81 


0. 


24 


-0 


.55 


-0. 


34 


258 


.24 


0.45 


257.79 


cs 2 


218. 


.26 


13 


.73 


2 


.04 


0. 


10 


-0 


.62 


-0. 


39 


233 


.12 


0.61 


232.51 


CH 2 


222. 


.18 


16 


.15 


1 


02 


0. 


52 


0. 


00 


-0 


06 


239 


.81 


0.59 


239.22 


CH 3 


204. 


53 


19. 


65 


1 


38 


0, 


49 





.00 


-0. 


05 


225 


.99 


-0.94 


226.93 


C 2 H 5 ( 2 A') 


175. 


.09 


11. 


.59 


-0 


.27 


0, 


20 





.00 


-0. 


06 


186 


.55 


-0.84 


187.39 


C3 H4 (cyclopropene) 


197. 


.38 


25. 


.54 


1 


.83 


0. 


05 


0. 


.00 


0. 


01 


224 


.82 


2.60 


222.22 


CH 2 CCH 2 


199. 


28 


24. 


.20 


1 


63 


0, 


19 





00 


0. 


02 


225 


.32 


1.59 


223.74 


sec-C3H7 


156. 


07 


15. 


35 





16 


0, 


.11 





00 


-0. 


06 


171 


.61 


-0.16 


171.77 


CgHg a 


187. 


.06 


25 


.50 


2 


29 




— 





00 




— 


214 


.85 


1.67 


213.18 


CN( 3 n cation) 


286. 


.12 


35 


.37 


2 


66 


0, 


89 





00 


-0. 


09 


324 


.95 


0.51 


324.44 


CN( 1 E+ cation) 


363. 


.28 


-31 


.31 


-12. 


94 


0, 


79 





.00 


0. 


02 


319 


.03 


0.11 


319.74 


CHO 


173. 


.77 


11 


.65 


0. 


50 


-0 


.01 


0. 


00 


-0. 


10 


185 


.81 


-2.14 


187.95 


CH 2 OH( 2 A) 


157. 


.75 


13. 


.94 


0. 


11 


0. 


19 


0. 


.00 


-0. 


08 


171 


.91 


-2.07 


173.98 


CH 3 0( 2 A') 


230. 


91 


18. 


.41 


-0 


,04 


0, 


07 





00 


-0. 


13 


249 


.22 


1.87 


247.34 


CH3OH 


231. 


.29 


21 


.46 





.84 


0. 


07 


0. 


00 


-0. 


11 


253 


.56 


2.55 


251.01 


CH 3 F 


279. 


.03 


14. 


.98 


-0. 


.55 


0, 


15 


0. 


00 


-0. 


09 


293 


.52 


3.85 


289.67 


CH 2 S 


190. 


95 


24 


.41 


1 


97 





09 





.00 


-0. 


26 


217 


16 


0.64 


216.52 


CH 2 SH 


161 


21 


11. 


50 


-0 


.62 


0, 


18 


0. 


.00 


-0, 


22 


172 


.04 


-2.35 


174.39 


CH3SH 


194. 


.71 


22. 


.63 


1 


58 


0, 


09 


0. 


.00 


-0. 


.25 


218 


.76 


0.83 


217.93 


CH3CI 


241. 


.14 


21 


17 


1 


07 


0, 


12 





00 


-0, 


30 


263 


.20 


2.92 


260.28 


CH3CHO 


207. 


52 


28. 


75 


1 


72 


0, 


22 





00 


-0, 


13 


238 


.08 


1.52 


236.56 


CH3OF 


243. 


.59 


20 


.99 


0. 


.43 


0, 


23 





00 


-0 


19 


265 


.04 


1.74 


263.30 


C 2 H4S(thiirane) 


183. 


16 


24. 


70 


1 


90 





15 





.00 


-0. 


27 


209 


.64 


0.62 


209.02 




294 


22 


ru 


00 


_2 




n 
u. 


81 
01 


n 
-u 


ns 

Uo 


n 

-u. 


12 






u.yo 


QfiQ zlR 


C4H 4 0(furan) a 


176. 


.04 


27 


.03 


2 


19 






0. 


.00 






205 


.27 


0.54 


204.72 


B 2 H 4 


191 


.08 


26 


.75 


2 


41 


0. 


55 


0. 


00 


0. 


00 


220 


.79 


0.07 


220.73 b 


NH 


294. 


73 


15. 


12 





85 


0, 


56 


-0 


.11 


-0. 


11 


311 


.03 


0.39 


310.64 


NH 2 


224. 


.17 


32. 


.83 


1 


99 


0, 


01 





00 


0. 


02 


259 


.02 


1.36 


257.66 


N 2 H 2 


206. 


.09 


16. 


32 





36 


-0 


.15 





.00 


0. 


01 


222 


.63 


1.26 


221.37 


N 2 H 3 


163. 


.56 


10 


.78 


-0 


39 


0, 


08 





.00 


-0. 


07 


173 


.95 


-1.66 


175.61 


HOF 


269. 


.88 


23. 


.22 


0. 


,87 


0, 


26 





00 


-0. 


19 


294 


.04 


0.44 


293.60 


SiHa^Ai) 


182. 


64 


25. 


81 


1 


29 


0, 


.57 


0. 


00 


0. 


13 


210 


.44 


-0.22 


210.66 


SiH 3 


173. 


72 


11. 


.57 


1 


.03 


0. 


34 


0. 


.00 


-0, 


13 


186 


.53 


-0.56 


187.09 


Si 2 H 2 


168. 


.88 


18. 


.88 


1 


,77 


-0. 


.07 


0. 


00 


-0 


33 


189 


.13 


0.22 


188.90 


Si 2 H 4 


153. 


71 


27. 


95 


3 


.66 


0. 


.45 





.00 


-0. 


06 


185 


.71 


-0.17 


185.88 


Si 2 H 6 


198. 


17 


23. 


84 


1 


.71 


0. 


42 





.00 


-0. 


20 


223 


.95 


1.44 


222.51 6 



a core correlation and scalar relativistic steps were skipped. 



W2 theory: Si 2 H 6 222.67 kcal/mol, B 2 H 4 220.75 kcal/mol 
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TABLE S— III. Components of Wl computed total atomization energies (kcal/mol) of the G2-1 



test set of molecules. 





SCF 


CCSD 


(T) 


Core 


spin-orbit 


Scalar rel. 


Wl 


ZPVE 


Wl 


Species 


limit 


limit 


limit 


corr. 


splitting 


effects 


TAE 




TAE0 


LiH 


34.30 


23.49 


0.00 


0.32 


0.00 


0.00 


58.10 


1.99 


56.11 


BeH 


50.44 


-0.52 


0.37 


0.51 


0.00 


-0.02 


50.78 


2.90 


47.87 


CH 


57.17 


25.96 


0.89 


0.14 


-0.05 


-0.04 


84.08 


3.97 


80.11 


CH 2 ( 3 Bi) 


154.73 


34.40 


0.97 


0.79 


-0.08 


-0.15 


190.67 


10.62 


180.05 


CHa^Ai) 


130.00 


48.99 


1.90 


0.37 


-0.08 


-0.08 


181.10 


10.24 


170.86 


CH 3 


243.39 


61.77 


1.91 


1.04 


-0.08 


-0.17 


307.86 


18.33 


289.53 


CH 4 


331.51 


85.05 


2.93 


1.21 


-0.08 


-0.19 


420.42 


27.56 


392.86 


NH 


48.59 


33.20 


1.21 


0.11 


0.00 


-0.07 


83.04 


4.58 


78.46 


NH 2 


116.80 


63.05 


2.56 


0.31 


0.00 


-0.15 


182.57 


11.71 


170.86 


NH 3 


203.27 


90.60 


3.96 


0.62 


0.00 


-0.25 


298.18 


21.17 


277.01 


OH 


69.28 


36.22 


1.75 


0.14 


-0.03 


-0.12 


107.25 


5.21 


102.04 


OH 2 


159.92 


69.62 


3.64 


0.36 


-0.22 


-0.26 


233.05 


13.15 


219.90 


FH 


99.88 


39.81 


2.27 


0.17 


-0.39 


-0.20 


141.56 


5.76 


135.80 


SiHa^Ai) 


113.49 


39.80 


0.67 


-0.11 


-0.43 


-0.22 


153.21 


7.17 


146.04 


SiH 2 ( 3 Bi) 


108.29 


24.91 


0.31 


-0.61 


-0.43 


-0.57 


131.90 


7.39 


124.51 


SiH 3 


182.66 


45.33 


0.53 


-0.45 


-0.43 


-0.62 


227.01 


13.07 


213.95 


SiH 4 


259.95 


64.61 


0.72 


-0.37 


-0.43 


-0.70 


323.78 


19.25 


304.53 


PH 2 


105.83 


47.19 


1.23 


0.15 


0.00 


-0.29 


154.10 


8.25 


145.85 


PH 3 


173.14 


66.89 


1.89 


0.21 


0.00 


-0.46 


241.67 


14.77 


226.90 


H 2 S 


133.57 


47.90 


2.18 


0.24 


-0.56 


-0.41 


182.92 


9.30 


173.62 


HC1 


80.80 


25.22 


1.47 


0.15 


-0.84 


-0.26 


106.54 


4.14 


102.39 


Li 2 


3.96 


20.25 


0.00 


0.17 


0.00 


0.00 


24.37 


0.48 


23.89 


LiF 


93.76 


41.75 


3.26 


1.02 


-0.39 


-0.22 


139.18 


1.27 


137.90 


C 2 H 2 


300.59 


94.53 


8.24 


2.40 


-0.17 


-0.27 


405.32 


16.68 


388.64 


C 2 H 4 


435.23 


119.62 


7.36 


2.30 


-0.17 


-0.33 


564.02 


31.47 


532.54 


C 2 Hg 


558.01 


146.90 


6.43 


2.34 


-0.17 


-0.39 


713.12 


45.97 


667.15 


CN 


86.30 


82.78 


10.29 


1.10 


-0.08 


-0.17 


180.22 


3.03 


177.18 


HCN 


205.43 


97.22 


9.04 


1.62 


-0.08 


-0.21 


313.03 


10.10 


302.93 


CO 


182.21 


68.72 


7.96 


0.91 


-0.31 


-0.14 


259.35 


3.11 


256.24 


HCO 


186.82 


83.20 


8.34 


1.13 


-0.31 


-0.26 


278.91 


8.00 


270.91 


0H 2 O 


265.28 


100.43 


7.89 


1.27 


-0.31 


-0.32 


374.25 


16.40 


357.85 


CH 3 OH 


376.88 


128.94 


6.98 


1.38 


-0.31 


-0.46 


513.41 


31.56 


481.85 


N 2 


120.96 


97.01 


9.34 


0.78 


0.00 


-0.11 


227.98 


3.45 


224.53 


H 2 NNH 2 


272.78 


156.60 


8.53 


1.10 


0.00 


-0.51 


438.50 


32.89 


405.61 


NO 


55.84 


86.80 


9.36 


0.44 


-0.05 


-0.16 


152.21 


2.78 


149.43 


o 2 


27.25 


83.68 


9.27 


0.25 


-0.45 


-0.15 


119.86 


2.29 


117.56 


HOOH 


142.56 


117.41 


9.12 


0.36 


-0.45 


-0.35 


268.66 


16.36 


252.30 


F 2 


-29.90 


61.02 


7.40 


-0.07 


-0.77 


0.03 


37.72 


1.48 


36.24 


co 2 


258.70 


116.21 


13.89 


1.69 


-0.53 


-0.45 


389.51 


7.23 


382.28 


Na 2 


-0.70 


17.25 


0.00 


0.58 


0.00 


0.02 


17.14 


0.22 


16.92 


Si 2 


44.86 


26.59 


4.30 


0.15 


-0.86 


-0.16 


74.88 


0.78 


74.11 


P2 


39.57 


66.81 


9.26 


0.63 


0.00 


-0.22 


116.04 


1.14 


114.90 


S 2 


49.49 


46.00 


7.45 


0.43 


-1.12 


-0.27 


101.97 


1.01 


100.97 


Cl 2 


26.55 


28.06 


4.40 


0.16 


-1.68 


-0.15 


57.34 


0.76 


56.58 


NaCl 


73.86 


24.02 


2.12 


-1.26 


-0.84 


-0.39 


97.51 


0.51 


97.00 


SiO 


118.16 


65.13 


8.28 


0.66 


-0.65 


-0.24 


191.34 


1.78 


189.56 


CS 


104.54 


56.57 


9.44 


0.72 


-0.64 


-0.15 


170.48 


1.84 


168.64 


SO 


53.01 


63.95 


8.41 


0.42 


-0.78 


-0.31 


124.69 


1.63 


123.06 


CIO 


9.58 


48.43 


6.26 


0.20 


-0.62 


-0.19 


63.65 


1.20 


62.46 


FC1 


15.17 


42.11 


5.30 


0.08 


-1.23 


-0.12 


61.31 


1.11 


60.21 


Si 2 H 6 


424.37 


109.43 


2.67 


-0.62 


-0.86 


-1.39 


533.61 


30.15 


503.46 


CH 3 C1 


303.82 


86.14 


5.29 


1.20 


-0.93 


-0.42 


395.10 


23.26 


371.84 


CH 3 SH 


358.26 


109.90 


6.07 


1.35 


-0.64 


-0.60 


474.35 


28.34 


446.01 


HOC1 


86.42 


72.76 


6.81 


0.28 


-1.06 


-0.28 


164.92 


8.14 


156.79 


so 2 


121.74 


120.94 


15.85 


0.83 


-1.01 


-0.71 


257.64 


4.35 


253.29 
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TABLE S-IV. Components of Wl computed total atomization energies (kcal/mol) of selected 
molecules from the G2-2 test set. 









U) 


Core 


spin-orbit 


Scalar rel. 


Wl 




Wl 


Species 


limit 


limit 


limit 


corr. 


splitting 


effects 


TAE 




1 Aii/U 


BF 3 


374.54 


87.43 


8.46 


1.76 


-1.18 


-0.68 


470.32 


7.65 


A O O C7fl. 

4oz.o7 


BCI3 


258.46 


56.11 


9.22 


1.60 


-2.55 


-0.76 


322.07 


4.64 


O "1 'I A AO. 

317.44 


CF4 


329.35 


136.15 


13.54 


1.04 


-1.63 


-0.82 


477.63 


10.42 


467.21™ 


OCS 


218.66 


101.05 


14.34 


1.61 


-0.87 


-0.53 


334.25 


5.71 


328.54 


CS2 


175.48 


87.14 


15.45 


1.52 


-1.20 


-0.63 


277.76 


4.29 


273.46 


N 2 


96.64 


154.02 


18.63 


1.18 


-0.22 


-0.41 


269.84 


6.93 


262.91 


C1NO 


55.24 


118.03 


17.09 


0.30 


-1.06 


-0.11 


189.49 


4.04 


"1 O FT A FT (1. 

185.45 


NF 3 


48.08 


141.59 


16.35 


-0.06 


-1.16 


-0.24 


204.57 


6.30 


198.27 


F2O 


-30.04 


109.50 


13.74 


-0.12 


-0.99 


-0.01 


92.07 


3.37 


88.70 


L>tl2r 2 


319.12 


110.47 


8.17 


1.06 


-0.85 


-0.52 


437.43 


20.19 


417.25 


CH3CN 


443.51 


157.83 


12.55 


2.86 


-0.17 


-0.41 


616.19 


27.92 


588.27 


HCOOH 


343.25 


144.99 


12.66 


1.59 


-0.53 


-0.58 


501.40 


20.85 


480.55 


{ > T T / 1 

O2H4U 


470.00 


166.92 


11.83 


2.41 


-0.39 


-0.55 


650.22 


35.36 


614.87 




437.15 


178.64 


17.09 


2.56 


-0.62 


-0.62 


634.20 


22.78 


611.41 


CH3OCH3 


597.05 


189.56 


10.63 


2.43 


-0.39 


-0.64 


798.63 


49.02 


749.60 


02H501 


533.26 


148.82 


9.13 


2.35 


-1.01 


-0.60 


691.95 


41.01 


a ten f\A 
OOU.94 


CH3COCH3 


736.37 


224.89 


15.49 


3.70 


-0.48 


-0.70 


979.26 


51.40 


927.86 


HS 


61.30 


25.22 


1.04 


0.12 


-0.02 


-0.19 


87.48 


3.76 


83.72 


CCH 


182.76 


73.53 


7.67 


2.08 


-0.17 


-0.27 


265.61 


8.68 


256.93 


C2H3 


339.57 


97.80 


6.70 


2.15 


-0.17 


-0.30 


445.75 


22.50 


423.25 


CH3CO 


424.39 


144.13 


11.92 


2.34 


-0.39 


-0.46 


581.93 


26.57 


555.36 


CH 2 OH 


294.61 


108.16 


6.54 


1.34 


-0.31 


-0.44 


409.90 


22.95 


386.95 


CH3CH2O 


524.99 


161.36 


9.33 


2.38 


-0.39 


-0.51 


697.16 


40.63 


656.53 


CH 3 S 


288.36 


88.31 


5.12 


1.28 


-0.64 


-0.39 


382.04 


21.83 


360.21 


C 2 H5 


472.09 


124.77 


5.54 


2.24 


-0.17 


-0.37 


604.09 


36.49 


567.60 


N0 2 


60.73 


146.46 


19.24 


0.73 


-0.45 


-0.40 


226.31 


5.43 


220.88 



(a) W2 values: BF 3 461.29 kcal/mol (C.W. Bauschlicher Jr., J. M. L. Martin, and P. R. 



Taylor, J. Phys. Chem. 103, 7715 (1999), using anharmonic zero-point energy of 7.89 
kcal/mol, Y. Pak and R. C. Woods, J. Chem. Phys. 106, 6424 (1997)); BC1 3 318.09 kcal/mol; 
CF 4 465.7 kcal/mol (estimated, see text); CS 2 274.35 kcal/mol; OCS 328.77 kcal/mol; C1NO 
185.46 kcal/mol; HS 83.70 kcal/mol; CCH 256.68 kcal/mol; F 2 88.49 kcal/mol; N 2 262.48 
kcal/mol; N0 2 220.42 kcal/mol. 
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TABLE S-V. Components (kcal/mol, at Wl level) of proton affinities of molecules in the G2/97 
test set 





SCF 


CCSD 


(T) 


Core 


Scalar rel. 


Final 


ZPVE 


Wl 


Thermal 


Final 


Species 


limit 


limit 


limit 


corr. 


effects 


Energy 




PA 


correction 


Wl PA 


NH 3 


216.00 


-3.12 


-1.02 


0.13 


-0.03 


211.94 


9.38 


202.57 


1.49 


204.06 


H 2 


175.57 


-3.13 


-0.92 


0.12 


-0.08 


171.56 


8.03 


163.53 


1.43 


164.97 


C 2 H 2 


163.47 


-4.61 


-1.22 


-0.27 


0.02 


157.39 


4.39 


153.00 


1.13 


154.13 


SiH 4 


157.46 


-0.04 


-0.47 


-0.28 


0.14 


156.82 


4.35 


152.47 


0.75 


153.22 


PH 3 


201.29 


-6.37 


-1.47 


-0.43 


-0.32 


192.70 


6.94 


185.76 


1.43 


187.19 


H 2 S 


176.84 


-1.72 


-0.84 


-0.07 


0.03 


174.24 


6.63 


167.62 


1.45 


169.07 


HC1 


137.91 


0.28 


-0.42 


-0.05 


0.03 


137.75 


5.17 


132.58 


1.17 


133.76 


H 2 


104.55 


1.72 


0.00 


0.00 


0.00 


106.27 


6.13 


100.14 


1.18 


101.32 



The first-order spin-orbit contribution is trivially zero in all cases. 
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